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We propose a new approach to the problem of rapidity gap survival (RGS) in the production of 
high-mass systems [H = dijet, heavy quarkonium, Higgs boson) in double-gap exclusive diffrac- 
tive pp scattering, pp ^ p + (gap) + H + (gap) + p. It is based on the idea that hard and soft 
interactions proceed over widely different time- and distance scales and are thus approximately 
independent. The high-mass system is produced in a hard scattering process with exchange of two 
gluons between the protons. Its amplitude is calculable in terms of the gluon generalized parton 
distributions (GPDs) in the protons, which can be measured in J/^p production in exclusive ep scat- 
tering. The hard scattering process is modified by soft spectator interactions, which we calculate 
in a model-independent way in terms of the pp elastic scattering amplitude. Contributions from 
inelastic intermediate states are suppressed. A simple geometric picture of the interplay of hard and 
soft interactions in diffraction is obtained. The onset of the black-disk limit in pp scattering at TeV 
energies strongly suppresses diffraction at small impact parameters and is the main factor in deter- 
mining the RGS probability. Correlations between hard and soft interactions {e.g. due to scattering 
from the long-range pion field of the proton, or due to possible short-range transverse correlations 
between partons) further decrease the RGS probability. We also investigate the dependence of the 
diffractive cross section on the transverse momenta of the final-state protons ("diffraction pattern"). 
By measuring this dependence one can perform detailed tests of the interplay of hard and soft in- 
teractions, and even extract information about the gluon GPD in the proton. Such studies appear 
to be feasible with the planned forward detectors at the LHC. 

PACS numbers: 12.38.-t, 13.85.-t, 13.85.Dz, 14.80.Bn 



I. INTRODUCTION 

Hard processes in high-energy pp scattering are impor- 
tant both as a laboratory for studying strong interaction 
dynamics and the parton structure of the proton, and as 
one of the main tools in the search for new heavy parti- 
cles. Of particular interest are so-called diffractive pro- 
cesses, in which the produced high-mass system (dijet, 
heavy particle) is separated from the projectile fragments 
by large rapidity gaps. Double-gap exclusive processes 
(i.e., without breakup of the protons) 

pp ^ p+ (gap) + H+ (gap) +p, (1) 

are considered as an option for the Higgs boson search 
at the LHC 0, S Hi S 0- Such processes have lower 
cross section than inclusive double-gap processes (with 
breakup of one or both protons), but offer better chances 
for detection, and for determining the mass of the pro- 
duced particle and possibly even its quantum numbers; 
see Ref. and references therein. Double-gap exclusive 
processes also appear to be an effective method for 
producing heavy quarkonia and investigating their prop- 
erties. 

From the point of view of strong interactions, double- 
gap exclusive events arise as the result of an interest- 
ing interplay of "hard" (involving momentum transfers 
much larger than the typical hadronic mass scale) and 
"soft" (momentum transfers of the order of the typical 



hadronic mass scale) interactions. The high-mass sys- 
tem is produced in a hard scattering process, involving 
the exchange of two gluons between the protons. The 
requirement of the absence of QCD radiation ensures the 
localization of this process in space and time. This alone, 
however, is not sufhcient to guarantee a diffractive event. 
One must also require that the soft interactions between 
the spectator systems do not lead to particle production. 
This results in a suppression of the cross section as com- 
pared to the hard scattering process alone, the so-called 
rapidity gap survival (RGS) probability. While not di- 
rectly observable, this quantity plays a central role in 
the discussion of hard diffractive processes and their use 
in new particle searches. 

Diffractive final states are most favorably produced in 
scattering at large impact parameters (peripheral scat- 
tering), where the chances for the spectator systems not 
to interact inelastically are large. Conversely, this means 
that the selection of diffractive events changes the effec- 
tive impact parameters as compared to inclusive events 
with the same hard scattering process. This effect is 
essential for understanding the physical mechanism of 
RGS. In this sense, RGS is a manifestation of a gen- 
eral quantum-mechanical phenomenon — the postselec- 
tion of certain initial-state configurations by conditions 
imposed on the final state 

The concept of RGS should also be viewed in the con- 
text of QCD factorization for the production of heavy 
particles in pp scattering . QCD factorization was for- 
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mally proved for inclusive scattering, pp — > H+X. A cru- 
cial element in this is the cancellation of initial-state and 
final-state QCD radiation. This cancellation becomes in- 
complete if additional conditions, such as rapidity gaps, 
are imposed on the hadronic final state. The introduc- 
tion of the RGS probability can be seen as an attempt to 
"restore" QCD factorization for diffractive processes at 
the phenomenological level, in a form analogous to the 
inclusive case. 

In this paper, we propose a new approach to RGS 
in double-gap exclusive hard diffractive processes. It is 
based on the idea that hard and soft interactions are 
approximately independent because they happen over 
widely different time- and distance scales. We imple- 
ment this idea in the framework of a partonic descrip- 
tion of the proton, along the lines of Gribov's parton 
picture of high-energy hadron-hadron scattering j^] . In 
the approximation where hard and soft interactions are 
considered to be completely independent, one can regard 
the hard scattering process as a "local operator" in par- 
tonic states, and RGS appears as the "renormalization" 
of this operator due to soft interactions. At the amplitude 
level, this leads to an absorption correction to the hard 
production process due to elastic rescattering; contribu- 
tions from inelastic intermediate states are suppressed 
because of the different character of the states accessible 
in hard and soft interactions. At the cross section level, 
we recover a simple "geometric" expression for the RGS 
probability, which was suggested on the basis of heuristic 
arguments in Refs. 9., 10, llj. In addition to providing a 
transparent physical picture of RGS, this expression can 
readily be evaluated in terms of two phenomenological 
ingredients, both of which can be probed in independent 
measurements: 

• The gluon generalized parton distribution (GPD) 
in the proton; more precisely its t-dependence 
("two-gluon formfactor" ) , whose Fourier transform 
describes the transverse spatial distribution of glu- 
ons. Information about it comes from measure- 
ments of hard exclusive processes in ep scatter- 
ing, in particular J/ijj photoproduction (HERA, 
FNAL). 

• The pp elastic scattering amplitude at high ener- 
gies; in particular its profile function in the impact 
parameter representation. It is known from fits to 
pp/pp total and elastic cross section data up to the 
Tevatron energy, and constrained by general argu- 
ments (black-disk limit, see below) 

The framework provided by our partonic approach to 
RGS allows us to take into account two basic facts about 
the dynamics of hard and soft interactions at high ener- 
gies, which turn out to have a decisive influence on the 
numerical value of the RGS probability. These are: 

• Small transverse radius of hard interactions. The 
radius of the transverse distribution of hard glu- 
ons in the proton is significantly smaller than the 



transverse radius of soft interactions in high-energy 
pp colHsions ("two-scale picture"). This basic fact 
explains many qualitative features of hard exclu- 
sive diffractive processes, such as the effective im- 
pact parameters in diffractive events, the order-of- 
magnitude of the RGS probability, and the pat- 
tern of the transverse momentum dependence of 
the cross section. 

• Black-disk limit in high-energy pp scattering. 
Parametrizations of the data as well as general the- 
oretical arguments indicate that the profile function 
of pp elastic scattering becomes "black" at small 
impact parameters at energies above the Tevatron 
energy, ^/s > 2 TeV. This circumstance makes 
the description of pp scattering at small impact 
parameters at LHC energies practically model- 
independent, and is essential for the stability of 
numerical predictions for the RGS. 

The evidence supporting these statements is described in 
detail in Ref. |9|, and summarized in Sections UTI and IlIII 
below. 

Our partonic approach also allows us to describe the 
effects of correlations between hard and soft interactions 
on the RGS probability. Such correlations can arise from 
various dynamical mechanisms, e.g. from the long-range 
pion field of the proton, or from possible short-range 
transverse correlations between hard partons, as sug- 
gested by the Tevatron CDF data on inclusive pp scatter- 
ing with multiple hard processes We find that the 
inclusion of such correlations decreases the RGS proba- 
bility compared to the independent interaction approx- 
imation. While these effects cannot be calculated in a 
completely model-independent way, they are important 
both for our general understanding of the mechanism of 
RGS, and for obtaining reliable numerical estimates of 
the RGS probability. These effects clearly merit further 
study. 

A unique feature of exclusive diffractive processes is 
that the interplay of hard and soft interactions can be 
studied experimentally, by measuring the dependence of 
the cross section on the transverse momenta of the final- 
state protons. The modification of the hard scattering 
amplitude by soft elastic rescattering can be viewed as 
an interference phenomenon, which gives rise to a dis- 
tinctive "diffraction pattern" in the final-state transverse 
momenta. By measuring this dependence in exclusive 
diffractive processes with relatively large cross section, 
such as dijet production, one can perform a variety of 
tests of the diffractive reaction mechanism, and extract 
information about the transverse radii of hard and soft 
interactions and their energy dependence. In Higgs pro- 
duction, measurements of the transverse momentum de- 
pendence would allow one to obtain additional informa- 
tion about the parity of the produced particle j^. Ex- 
perimentally, such studies appear to be feasible with the 
planned forward detectors at the LHC 1^ and the Teva- 
tron Q. 
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Detailed studies of RGS in diffractive pp scattering 
were done in a model of soft interactions based on eikon- 
alized Pomeron exchange 0,1^. We show here that, in 
the approximation where hard and soft interactions are 
considered to be independent, the RGS probability un- 
ambiguously follows from QCD and can be calculated 
in a model-independent way. Nevertheless, our numeri- 
cal results for the RGS probability in this approximation 
turn out to be of the same order of magnitude as those 
reported in Refs. 0,01) which can be attributed to the 
fact that the Pomeron parametrization of the pp elastic 
amplitude reproduces the approach to the BDL at small 
impact parameters. The agreement between our numer- 
ical results and theirs at the quantitative level is some- 
what accidental, being due to the fact that in the calcu- 
lation of Refs. Qilal the effect from inelastic intermediate 
states (which in our approach are seen to be strongly sup- 
pressed because of the small overlap of states accessible 
via hard and soft interactions) is partly compensated by 
the choice of a larger value of the t-slope of the gluon 
GPD; see Section[V]for details. Finally, our partonic ap- 
proach can naturally be extended to include correlations 
between hard and soft interactions, which have a poten- 
tially large numerical effect on the RGS probability. 

This paper is organized as follows. In Section ^ we 
review the information about the transverse structure 
of soft interactions from pp elastic scattering, and the 
approach to the black-disk limit at central impact pa- 
rameters. In Section irm we discuss the properties of the 
proton's gluon GPD at small x and summarize our knowl- 
edge of the transverse spatial distribution of hard partons 
in the proton. Section IIVI describes the basic framework 
of our approach to RGS. We outline the properties of 
the hard scattering amplitude, describe the theoretical 
formulation of the independence of hard and soft inter- 
actions, and obtain a master expression for the diffractive 
amplitude combining hard and soft interactions. We then 
explain the suppression of inelastic intermediate states, 
and evaluate the diffractive amplitude in terms of the 
gluon GPD and the pp elastic amplitude. In Section ]V\ 
we use our result for the amplitude in the independent 
interaction approximation to calculate the RGS proba- 
bility. We recover a simple geometric expression for the 
RGS probability and discuss the effective impact param- 
eters in exclusive diffraction. We then evaluate the RGS 
probability numerically, estimate the uncertainty of the 
numerical predictions due to the phenomenological input, 
and emphasize the crucial role of the black-disk limit in 
stabilizing the numerical predictions. We also comment 
on the results for the RGS obtained within the eikonal- 
ized Pomeron model for soft interactions 4] from the per- 
spective of our approach. In Section IVII we discuss var- 
ious effects beyond the approximation of independence 
hard and soft interactions in exclusive diffraction. We 
first point out that hard screening corrections may re- 
duce the diffractive cross section beyond the RGS proba- 
bility due to soft interactions. We then discuss the effect 
of correlations between hard and soft interactions on the 



RGS probability, considering two specific mechanisms — 
diffractive scattering from the long-range pion field, and 
short-range transverse correlations between partons. In 
Section rVlII we work out the dependence of the exclusive 
diffractive cross section on the final proton transverse 
momenta. We discuss which experimentally observable 
features of this dependence furnish useful tests of the 
diffractive reaction mechanism, and how one can extract 
information about the gluon GPD. In Section IVIIII we 
summarize our results. We comment on the implications 
for the Higgs boson search, and on the experimental fea- 
sibility of measuring the transverse momentum depen- 
dence of exclusive diffraction with the planned forward 
detectors at the LHC. 



II. BLACK-DISK LIMIT IN pp ELASTIC 
SCATTERING 

Information on the transverse radius of strong interac- 
tions at high energies comes mostly from measurements 
of the ^-dependence of the differential cross section for 
pp and pp elastic scattering. Combining these data with 
those on the pp/pp total cross section, and implement- 
ing theoretical constraints following from the unitarity of 
the S'-matrix, one can reconstruct the complex pp elastic 
scattering amplitude, Tci{s,t); see e.g. Refs. J^,16„Jj|. 
At high energies, s ^ \t\, angular momentum conser- 
vation in the CM frame implies that the scattering am- 
plitude is effectively diagonal in the impact parameter of 
the colliding pp system. It is convenient to represent it as 
a Fourier integral over a transverse coordinate variable, 
b, 



Tci{s, t 



) = ^ y'd^6e-(^-'')r(.,6), (2) 

where F is the (dimensionless) profile function. One can 
then express the elastic, total, and inelastic (total minus 
elastic) pp cross sections in terms of the profile function 

as 



a'tot(s) 

O-incl(s) 



d^b X 



2ReF(s,6), 

[i-|i-r(.,b)|2] 



(3) 



The functions on the R.H.S. describe the distribution of 
the respective cross sections over pp impact parameters, 
6 = |b| >Qu\. In particular, we note that the combination 

r(.,6)p (4) 
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can be interpreted as the probability for "no inelastic in- 
teraction" in a pp collision at impact parameter 6; this 
combination plays an important role in our calculation of 
the RGS probability (see Section M below) |6lli |. A mea- 
sure of the transverse size of the proton is the logarithmic 
t-slope of the elastic pp cross section at i = 0, 



dt 



d(7c\/dt (t) 



dac\/dt{Q) 



(5) 
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At high energies, where the elastic amphtude is predom- 
inantly imaginary, and T is real, B is equal to half the 
average squared impact parameter in the total pp cross 
section, 



B 



_ Ij (fbb^2Rer{s,b) 
2 J d26 2Rer(s,6) 



(6) 



which may be associated with the transverse area of the 
individual protons. The data show that the slope in- 
creases with the CM energy as 



B{s) = B{sq) + 2a' ln(s/so), 



(7) 



where a' « 0.25 GeV^^. In the Pomeron exchange 
parametrization of the pp elastic amplitude this constant 
is identified with the slope of the Pomeron trajectory. 

In Gribov's parton picture of high-energy hadron- 
hadron interactions Q , the transverse size of the proton 
in pp elastic scattering can be directly associated with 
the average transverse radius squared of the distribution 
of soft partons mediating the soft interactions, 



B = (p')soft. 



(8) 



Here and in the following, we use p = \p\ to denote the 
transverse distance of partons from the center of the pro- 
ton, and 6 = |b| for the impact parameter of the pp col- 
lision. The growth of the proton's transverse size with 
energy is explained as the result of random transverse dis- 
placements in the successive decays generating the distri- 
bution of soft partons (Gribov diffusion). Below we shall 
compare this distribution of soft partons to the distribu- 
tion of hard partons probed in hard exclusive processes 
(see Section llll|) . 

Parametrizations of the available data indicate that at 
energies above the Tevatron energy, > ^/sTcvation = 
2 TeV, the profile function at small impact parameters 
approaches 



r(s,6) 



1 



for b < 6o(s). 



(9) 



This corresponds to unit probability for inelastic scatter- 
ing for impact parameters b < 6o(s), cf. Eqs. © and Q), 
similar to the scattering of a pointlike object from a black 
disk of radius and is referred to as the black-disk limit 

(BDL) UMm 

The approach to the BDL in central pp scattering at 
high energies is a general prediction of QCD, independent 
of detailed assumptions about the dynamics. Studies of 
the interaction of small-size color dipoles with hadrons, 
based on QCD factorization in the leading log approx- 
imation, show that the BDL is attained at high energies 
as a result of the growth of the gluon density at small x 
due to DGLAP evolution jH. This result can be used to 
estimate the interaction of leading projectile partons with 
the small-x gluons in the target in pp scattering; one finds 
that there is no chance for the projectile wavefunction to 
remain coherent in small impact parameter scattering at 
TeV energies lllj . Similar reasoning allows one to 
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FIG. 1; The probability distribution for no inelastic inter- 
action, |1 — r(s,f))|^, Eq. as a function of fe = at 
the LHC energy (y^ = 14 TeV), as computed with different 
parametrizations of the pp elastic scattering amplitude. Solid 
line: Parametrization of Islam et al. \vi\ ("diffractive part" 
only). Dashed line: Exponential parametrization, Eq. (I12L 
with Fo = 1 (BDL), cf. Eq. lO, and B = 21.8 GeV"^ 



predict t he g rowth of the size of the black region, bo, 
with s ^U^. As a by-product, these arguments explain 
why the observed coefhcient in the Froissart formula for 
the total cross sections is significantly smaller than that 
derived from the general principles of analyticity of the 
amplitude in momentum transfer and unitarity of the S- 
matrix J^J. We note that the need for the approach to 
the BDL in high-energy scattering at central impact pa- 
rameters was understood already in the pre-QCD period 
within the Pomeron calculus, where it was noted that 
this phenomenon resolves the apparent contradiction be- 
tween the formulae of the triple-Pomeron limit and the 
unitarity of the ^-matrix, especially in models where the 
Pomeron intercept, ap{0), exceeds unity pol |. 

For our studies of diffractive pp scattering it will be 
useful to have a simple analytic parametrization of the 
pp elastic amplitude at the LHC energy, which incorpo- 
rates the approach to the BDL at small impact parame- 
ters. The t-dependence of the pp elastic scattering cross 
section for \t\ ^ 1 GeV^ over the measured energy range 
is reasonably described by an exponential shape. 



(X exp [B{.s) t] , 



(10) 



where B(s) represents an effective slope, to be distin- 
guished from the "exact" slope at < = 0, Eq. JSJ. A 
parametrization of the pp elastic amplitude which repro- 
duces this dependence is 



corresponding to 

r{s,b) 



o-tot(s) exp 

OTT 



Fo(s) exp 



B{s)t 



2B{s) 



(11) 



(12) 
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with 

ro(.) . r(.,6^o) ^ £1^. (13) 

Equation pi|) takes into account that the amplitude at 
high energies is predominantly imaginary, and satisfies 
the optical theorem for the total cross section, (Ttot(s) = 
(Stt/s) ImToi(s, i = 0). We may now incorporate the 
constraint of the BDL at small impact parameters by 
replacing 

To -1. (14) 

The value of B we determine by comparing the profile 
function H12I) with phcnomcnological parametrizations of 
the data, extrapolated to the LHC energy, which gives 

B w 20GeV"^ (\/i=14TeV). (15) 

In particular, with B ~ 21.8 GeV^^ we obtain excellent 
agreement with the Regge parametrization of Ref. 0. 
Figure^shows the probability for no inelastic interaction, 
|1 — r(s,6)p, Eq. computed with the phenomeno- 
logical parametrization of Ref. and our exponential 
parametrization incorporating the BDL, Eqs. 1(12(1 and 
(114(1 . One sees that the simple exponential parametriza- 
tion is a reasonable overall approximation to the phe- 
nomcnological parametrization over the 6-range shown 
in Fig. CI 

III. TRANSVERSE SPATIAL DISTRIBUTION 
OF GLUONS 

Information about the transverse structure of hard 
interactions comes from studies of hard exclusive pro- 
cesses in ep scattering, such as meson electroproduction 
or virtual Compton scattering. Such processes probe the 
GPDs in the proton, whose Fourier transform with re- 
spect to the transverse momentum transfer to the proton 
describes the spatial distribution of quarks and gluons in 
the transverse plane; see Refs. for a review. In 

this section we summarize what is known about the gluon 
GPD at small x from theoretical considerations, and from 
measurements of J/ip photoproduction and other pro- 
cesses at HERA and in fixed-target experiments. 

The gluon GPD can be formally defined as the transi- 
tion matrix element of the twist ~2 QGD gluon operator 
between proton states of different momenta, p and p' . 
Physically, it describes the amplitude for a fast-moving 
proton to "emit" and "absorb" a gluon with given lon- 
gitudinal momenta, with transverse momenta (virtual- 
ities) integrated over up to some hard scale, Q^, and 
a certain invariant momentum transfer to the proton, 
t = {p' —p)^. The choice of longitudinal momentum vari- 
ables is a matter of convention. Instead of the initial and 
final gluon momentum fractions (with respect to the ini- 
tial proton momentum), x and x\ we use as independent 
variables the initial gluon momentum fraction, x, and the 



fractional lon gitu dinal momentum transfer to the proton 
("skewness") [63 

f = x~x', (16) 
and denote the gluon GPD by 

i/g(x,^,t;Q2). (17) 

In the limit of zero momentum transfer, the gluon GPD 
reduces to the usual gluon momentum density in the pro- 
ton m 

Hg{x,^ = Q,t^O;Q^) = xG{x,Q^). (18) 

For discussing the ^-dependence of the gluon GPD it is 
convenient to write it in the form 

Hg{x,tt-Q^) = Hg{x,^,t^O;Q^) Fg{x,^,t;Q^), 

(19) 

where the function Fg is known as the "two-gluon form- 
factor" of the proton and satisfies Fg{t — 0) ~ 1. Note 
that the two-gluon formfactor still depends on x and ^, 
i.e., Eq. ((19(1 does not imply naive factorization of the 
i-dependence from that on the partonic variables. 

The dependence of the gluon GPD on the QCD scale, 
Q^, is governed by the QCD evolution equations. In ap- 
plications to production of fixed-mass systems in high- 
energy ep or pp collisions with <C s (such as Higgs 
boson production at the LHC) we shall be interested in 
the gluon GPD in the region where 

x,x' < 1, (20) 

while at the same time is much larger than the typ- 
ical hadronic mass scale, ^ IGeV^. In this region 
the gluon GPD can be calculated by applying QCD evo- 
lution to a "primordial" distribution at a low scale, Qq, 
in which one neglects the skewness, ^ = 0, or x = a:' (di- 
agonal approximation) 2,^,2^- QCD evolution degrades 
the individual gluon momentum fractions with increasing 
Q^, while their difference, ^ = x — x' , remains fixed by 
kinematics, being equal to the longitudinal momentum 
transfer to the proton; as a result, the primordial GPD 
at the low scale is sampled in the region jcc — a;'| <^ x, x' , 
where the diagonal approximation is justified. In the di- 
agonal approximation, the GPD at t = is completely 
determined by the usual gluon density, of. Eq. ((TH|l . leav- 
ing only the t-dependence (two-gluon formfactor) and 
its correlation with x up to modeling. This approxima- 
tion makes for a great simplification in applying GPDs 
at small x, and has been used extensively in the analysis 
of exclusive electroproduction processes; see Ref. 9] for 
a review. 

For X 1, the two-gluon formfactor permits a sim- 
ple interpretation in terms of a spatial distribution of 
gluons in the proton. For f <C 1, the invariant momen- 
tum transfer is dominated by the transverse momentum 
transfer between the proton states, 

t « -Ai, = pl-pi. (21) 
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The two-gluon formfactor can be represented as a Fourier 
integral over a transverse coordinate variable, p. 



P e 



X Fg{x,i,p-Q^). (22) 

For economy of notation we use the same symbol for the 
two-gluon formfactor and its Fourier transform, distin- 
guishing the two functions by their arguments. The p- 
dependent function describes the spatial distribution of 
gluons in the proton in the transverse plane; see Ref. 
for a review. For ^ = (i.e., a;' = x) it is positive definite 
and can be interpreted probabilistically as the gluon den- 
sity at transverse position p 26] ; for f 7^ it describes 
the non-diagonal transition matrix element of the gluon 
density j27i ] . A measure of the gluonic transverse size of 
the nucleon for given x and is the average of p^, which 
is proportional to the t-slope of the two-gluon formfactor 
at t = 0, 



d^pp" Fg{x,i = 0,p;Q^) (23) 



4^ 

dt 



X, e = 0, i = 0; Q2 



(24) 



The two-gluon formfactor of the nucleon, and hence 
the transverse spatial distribution of gluons, can directly 
be extracted from the i-dependence of the differential 
cross section for hard exclusive vector meson production 
processes probing the gluon GPD. QCD factorization im- 
plies that the i-dependence of the cross section resides 
solely in the gluon GPD, 



day 

dt J 



cx 



F^{x,^,t; 



(25) 



up to small higher-twist corrections due to the finite size 
of the produced vector meson |^. In particular, the t- 
slope at t = is proportional to the proton's average 
gluonic transverse size. 



d 

dt 



da/dt{t) 



da/dt (0) 



(26) 



t=o 



A crucial test of the applicability of QCD factorization is 
provided by the observed convergence of the t-slopes of 
various gluon-dominated vector meson production pro- 
cesses {J/ip,p,4>) at large Q^; see Ref. [3 for a detailed 
discussion. 

A particularly clean probe of the two-gluon formfac- 
tor is J/ip photoproduction, the i-dependence of which 
in experi 

(x ~ 10-2 - 10-4) liSlHil 

periment {{x) ^ 5 x IQ-^) [23], and a number of fixed- 
target experiments at lower energies {x ~ 10^^); see 
Refs. 0, for a review of the data. This process 
probes the two-gluon formfactor at an effective scale 
Q2 ^ SGeV^. Analysis of the data, combined with the- 
oretical investigations, has produced a rather detailed 



has been measured in experiments at the HERA collider 

ij, the FNAL fixed-target ex- 
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FIG. 2: The "two-scale picture" of the transverse structure 
of the proton in high-energy collisions. 



picture of the gluonic transverse size of the nucleon 
and its x-dependence j^]- For x ^ 0.1 — 0.3, the glu- 
onic transverse size suggested by the fixed-target data is 
~ 0.25 fm^, close to 2/3 times the proton's axial 
charge radius, (r^)^. Between x ~ IQ-^ and x ~ 10-^, 
(p2)g increases by ~ 30%. This can be explained by the 
contribution of the nucleon's pion cloud to the gluon den- 
sity at large transverse distances, p ^ 1/{2M.^), which is 
dynamically suppressed for x > M^/Mj^j and reaches its 
full strength for x <C Mt^/Mn 33]. Finally, over the 
HERA range, x ~ IQ-^ — 10^*, the gluonic transverse 
size exhibits a logarithmic growth with 1/x, 

{p')9 = (p')g(xo) + 4a; In^ (x < xo « 10-2), 

(27) 

with a rate, a^, considerably smaller than that governing 
the growth of the proton's transverse size in pp elastic 
scattering, which is dominated by soft interactions, 



a'g < a'. 



(28) 



A recent analysis of the HI data finds a'g = 0.164 ± 
0.028(stat)±0.030(syst) GeV^^ for J/V^ photoproduction 
and 0.019 ± 0.139(stat) ± 0.076(syst) GeV-^ for electro- 
production [23; an analysis of ZEUS electroproduction 
data quotes a'g = 0.07±0.05(stat)t;5:o4(syst) GeV^^ ;3(|, 
significantly smaller than the soft value a' = 0.25 GeV^^. 
The smaller rate of growth of the nucleon's size in hard 
interactions can qualitatively be explained by the sup- 
pression of Gribov diffusion in the decay of hard (highly 
virtual) partons as compared to soft partons. 

A crucial observation is that the transverse area occu- 
pied by partons with x ^ IQ-^ is much smaller than the 
transverse area associated with the proton in soft inter- 
actions (see Fig.|21l, 



{p'),{x>lO-') « 



(P^)soft, 



or 



2Bg < B. 



(29) 



(30) 



In high-energy pp collisions with hard partonic processes 
one is thus dealing with a two-scale picture of the trans- 
verse structure of the proton. Moreover, when consider- 
ing the production of a heavy particle with fixed mass. 
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Dipole, = 1 GeV^ 
- \ Exp., B = 3.24 GeV'- 




corresponding to 



exp [- 



pV(2i?,) 



2ttB„ 



(34) 



The relation between the parameters of the dipole and ex- 
ponential parametrization which would follow from iden- 
tifying = 4:dFg/dt{t = 0) is Bg = A/ml. Better 
overall agreement between the squared formfactors for 
\t\ < 1 GeV^ is obtained for somewhat smaller values of 
Bg. Matching the squared formfactors at |i| = 0.5 GeV^ 



we obtain 



-t [GeV"] 

FIG. 3: Comparison of the dipole (solid line) and exponential 
(dashed line) parametrizations of the two-gluon formfactor, 
with the parameters related by Eq. (1351 . Shown is the squared 



B„ 



3.24 



(35) 



see Figure |21 It was shown in Ref. that both the 
dipole with m^ = 1.1 GeV^ and the exponential with 



two-gluon formfactor, Fg{t), for both parametrizations, as ^ 3.0 GeV^ given by Eq. ^ describe well the t- 



corresponds to the f-dependence of the cross section for J/t/j 
photoproduction (details see text) 



TO_f/, in a partonic process with xi.2 ^ mn j \fs, the "soft" 
area of the proton increases with s faster than the "hard" 
area (which changes as a result of the decrease of x) , be- 
cause o! > a'g, cf. Eq. (|28|l . Thus, the difference of the 
two areas becomes even more pronounced with increasing 
energy. 

For our studies of hard processes in diffractive pp scat- 
tering we require a parametrization of the i-dependence 
of the two-gluon formfactor, viz. the shape of the trans- 
verse spatial distribution of gluons. The a:-values probed 
in Higgs production at central rapidities are x ^ 10~^ 
at the LHC energy. Taking into account the effect of 
DGLAP evolution, even larger values of x are probed 
when parametrizing the two-gluon formfactor at the J/4' 
production scale, ~ 3 GeV^ (for a general discussion 
of the effect of DGLAP evolution on the transverse spa- 
tial distribution of gluons, see Ref. 0|). We thus need 
to look at the J/ip photoproduction data at a; ^ 10^^, 
which are probed in fixed-target experiments. 

Theoretical arguments suggest that the two-gluon 
formfactor at a; > 10~^ should be close to the axial form- 
factor, which is well described by a dipole form (we omit 
all arguments except t), 



1 



{i-t/mir- 



(31) 



1 GeV^ • The corresponding transverse 



with to; 

spatial distribution of gluons is given by 

2 

rr 



a 
27r 



(^) K,{mgp). (32) 



We also consider an exponential parametrization of the 
two-gluon formfactor. 



Fg{t) = exp(Sgt/2), 



(33) 



dependence of the data from the FNAL E401/E458 ex- 
periment at {E^) = 100 GeV in which the recoiling pro- 
ton was detected [2^ . We also note that this value of Bg 
is consistent with what one obtains from the extrapola- 
tion of the HERA data towards larger x, using Eq. ()27|) 
with the measured a'g. We shall use the dipole, Eqs. (|31|l . 
with TOg = 1 GeV^ and the exponential, Eq. (|33|l , with 
Bg — 3.24 GeV~^, as our standard parametrizations for 
calculations in the kinematics of Higgs production at the 
LHC below; comparison between the two will allow us 
to estimate the uncertainty of our numerical predictions 
with respect to the shape of the two-gluon formfactor. 



IV. THEORY OF RAPIDITY GAP SURVIVAL 

We now outline the basic steps in the calculation of the 
amplitude of double-gap exclusive diffractive processes 
and develop the physical picture of RGS. The un- 
derlying idea of our approach is that hard and soft in- 
teractions are approximately independent, because they 
happen over widely different distance- and time scales. 



A. Hard scattering process 

In the first step, one calculates the amplitude for 
double-gap diffractive production of the high-mass sys- 
tem due to hard interactions. For definiteness, we shall 
refer in the following to Higgs boson production, keeping 
in mind that the discussed mechanism applies to pro- 
duction of other high-mass states as well (dijets, heavy 
quarkonia, etc.). According to electroweak theory, the 
Higgs boson is produced predominantly through its cou- 
pling to gluons via a quark loop; for a review and ref- 
erences see Ref. [s^]. In contrast to inclusive produc- 
tion, the amplitude for double-gap diffractive produc- 
tion is in the lowest order in the QCD running coupling 
constant, a^, given by the exchange of two gluons with 
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vacuum quantum numbers in the ^-channel (see Fig.^. 
The Higgs boson is radiated from one of the gluon hues. 
The role of the second exchanged gluon is to neutralize 
the color charge in order to avoid gluon bremsstrahlung. 
However, global color neutrality alone is not sufhcient. 
To suppress radiation, one must require that color be 
screened locally in space-time. Conversely, this means 
that the selection of a diffractive process, without accom- 
panying radiation, guarantees some degree of localization 
of the exchanged system. 

Operationally, the localization of the exchanged two- 
gluon system is ensured by Sudakov formfactors, which 
suppress configurations with low virtualities prone to 
emit gluon bremsstrahlung. The actual calculation of the 
hard scattering amplitude including Sudakov suppression 
is a challenging problem, which was addressed in various 
approximations in Refs. |35l l3^ . Fortunately, for our 
purposes we do not need to solve this problem at a fully 
quantitative level, as only a few qualitative aspects of the 
hard scattering process turn out to be essential for the 
physics of RGS. 

To discuss the hard scattering process, it is natural to 
perform a Sudakov decomposition of the four-momenta, 
using the initial proton momenta, pi and p2, as basis 
vectors, with 2{pip2) = s (we neglect the proton mass). 
As the transverse momenta of the final-state protons are 
small compared to the Higgs mass, we can expand the 
final proton four-momenta as 



Pi 

P2 



(1 
(1 



6)Pi + P'i±, 

6)P2 + P'2±: 



(36) 



where {p[j^,pi) = {p[^,p2) = etc., and ^1^2 parametrize 
the longitudinal momentum loss [cf. Eq. H16|) and the 
footnote before it]. 



rriH 



(37) 



where y is the rapidity of the produced Higgs boson. As- 
suming a Higgs mass of the order of 100 — 200 GeV, the 
typical values of ^1,2 are of the order of 10^^ for produc- 
tion at central rapidities at the LHC {^/s — 14TeV). 

Consider now the two-gluon exchange process of Fig.^ 
as a Feynman diagram, in which the upper and lower 
blobs denote the gluon-proton scattering amplitudes, to 
be specified in more detail below. The four-momenta of 
the gluons coupling to the Higgs we parametrize as 



fci 

fc2 



XlPl 

x[pi 



X2P2 
X2P2 



-p'l^, 



(38) 
(39) 



The four-momentum of the screening gluon then follows 
from four-momentum conservation. 



k — x'lPi 



where 



X'2P2 



£.1 - xi, 

6 + X2. 



(40) 

(41) 
(42) 




FIG. 4: The hard scattering process in double-gap exclu- 
sive diffractive Higgs boson production Q. Two gluons are 
exchanged between the protons. The gluon-Higgs coupling 
is indicated as a local vertex. The upper and lower blobs 
denote the gluon-proton scattering amphtude, which can be 
calculated in terms of the gluon GPD in the proton. 



We want to identify the dominant region of integration 
in the loop integral. First, analogy with inclusive pro- 
duction of heavy particles at central rapidities suggests 
that the momentum fractions of the gluons producing the 
Higgs boson (with respect to their parent protons) are 
practically the same as given by the naive parton model 
estimate, in which one neglects the transverse momenta 
and virtualities of the annihilating gluons. 



™H ±y 

xi,2 ~ ■ 
\ s 



(43) 



That is, the momentum fractions of the annihilating glu- 
ons are equal to the protons' fractional longitudinal mo- 
mentum loss, Xi^2 ~ £.1,2- Second, in the case of double- 
gap diffractive production the Sudakov formfactor associ- 
ated with the Higgs boson vertex, which accounts for the 
absence of gluon bremsstrahlung, restricts the (spacelike) 
virtualities of the annihilating gluons and their transverse 
momenta to values of the order of some "intermediate" 
hard scale. 



h,i2, fi-J, 



)2 

C int 5 



with 



m^/4 > Qf„t > a: 



^QCD- 



(44) 



(45) 



This was demonstrated explicitly in Ref. |35|. where 
the distribution of transverse momenta in the loop in- 
tegral was studied in a model which included the LO 
Sudakov formfactor associated with the ggH vertex; see 
also Ref. js^ [b^I . Expressing now the virtualities of the 
annihilating gluons, k"^ 2^ ™^ terms of the decompositions 
(|38|1 and (|39|l . neglecting the proton transverse momenta 
relative to k±, we find that Eq. (|44|l implies 



'^2,1 



X2,lS 



Qfnt 
X2,lS 



Qint 



< 2^1,2, 



(46) 
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i.e., the energy and longitudinal momentum fraction of 
the screening gluon are substantially smaller than those 
of the annihilating gluons. The screening gluon does not 
"belong to" any of the two protons; its momentum is 
predominantly transverse, and it has spacelike virtuality, 

-e « -kl ^ gf„t- (47) 

In the annihilating gluons, on the other hand, longitudi- 
nal and transverse momenta contribute in equal amounts 
to the virtuality, 

~ xi,2x'2,is - kl - Ql,. (48) 

To summarize, the hard scattering process takes the form 
of the exchange of two gluons with comparable virtuali- 
ties ~ Qfnt' transverse momenta Qint, between the 
protons. Of the two gluons, one carries substantial lon- 
gitudinal momentum fraction of the proton, ~ mH/\/s, 
and annihilates with the corresponding other to make the 
Higgs, the other gluon represents a "Coulomb-like" ex- 
change with small momentum fraction Q^^^/ {mH\/s)- 
The important point about the two-gluon exchange 
process is the appearance of the intermediate hard scale, 
Qf^i, governing the virtualities and transverse momenta 
of the exchanged gluons. This allows us to make a cru- 
cial simplification in the description of the gluon-proton 
scattering amplitudes, which we have not yet specified 
so far. Namely, we argue that, in a partonic description 
of the proton, the gluon-proton scattering amplitude is 
dominated by the two gluons coupling to the same parton 
■ This approximation is analogous to the assumption 
of dominance of the "handbag graph" in virtual Comp- 
ton scattering at large photon virtuality, Q^, which is well 
established and forms the basis of QCD factorization for 
this process |3^ ^3 ^> 03 • In this approximation, 
the gluon-proton scattering amplitude, which is predom- 
inantly imaginary at high energies, can be calculated in 
terms of the generalized parton distributions (GPDs) — 
here, predominantly, gluon distributions — in the pro- 
tons. 

We do not attempt to calculate the absolute normal- 
ization of the amplitude for double-gap hard diffractive 
production through two-gluon exchange in terms of the 
GPDs; doing so would require a substantially more ac- 
curate evaluation of the two-gluon exchange graph than 
the qualitative estimates presented above. Fortunately, 
for the theory of RGS, the only information we require (in 
addition to the qualitative properties of the hard process 
derived above) is the dependence of the double-gap hard 
diffractive amplitude on the transverse momentum trans- 
fers to the protons, pi_L and P2±- For sufficiently large 
scales QfntJ this dependence is described by the GPDs, 
even if the gluon momentum fractions in the hard ampli- 
tude and the virtuality of the exchanges are subject to 
integration over a certain range, and determined only in 
order-of-magnitude. Thus, we can state that the pij,- 
and p2i_ dependence of the double-gap hard diffractive 



amplitude is proportional to 

Thard (X Hg{xi,^i,ti;Q^) Hg{x2,^2,t2;Q^), (49) 

where 

h = p[\ < 0, t2 - p'i^ < (50) 

are the invariant momentum transfers to the proton. 
Here the longitudinal momentum transfers, ^1^2, are kine- 
matically fixed by Eq. H37|l . while parton momentum frac- 
tions, Xi^2-, are determined by Eq. (|43|l with accuracy 
given by Eq. (|45|l . The resolution scale, Q^, at which 
the GPD needs to be taken here, is parametrically of the 
order Qf^^ [cf. Eq. but numerically substantially 

larger, 

^ several times Qf^f (51) 

This follows from the fact that, by convention, de- 
termines the upper limit of the transverse momenta in 
the parton distribution, while Qf^^ is a measure of the 
dominant (average) values in the distribution, which for 
a distribution is significantly lower than the upper 

limit. 

To conclude this discussion, some comments concern- 
ing QCD evolution and the modeling of the GPDs are 
in order. First, in Higgs boson production at the LHC 
we are dealing with the gluon GDP at a;, ^ ~ 10~^ and 
x' <C X, where it can legitimately be calculated by apply- 
ing QCD evolution to a diagonal GPD at a lower scale 
(cf. the discussion in Section ITTT|l . Second, at LHC en- 
ergies the typical momentum fraction of the screening 
gluon reaches values x[ 2 ~ 10~^ for ~ few GeV^, cf. 
Eq. . For such low values of x the use of DGLAP 
evolution in principle requires justification. However, as 
explained in detail in Ref. for such values of x the 
kinematic conditions still restrict the actual number of 
radiated gluons to a few, so that NLO DGLAP and re- 
summcd BFKL evolution give similar results, see Ref. 
for a review. In this sense, the use of GPDs generated by 
DGLAP evolution seems to be appropriate. Third, one 
may wonder about the breakdown of the pQCD calcula- 
tion of the hard reaction amplitude due to the growth of 
the gluon density at small x within the DGLAP approx- 
imation — the BDL, see Ref. Q and references therein. 
This should not be a problem in the present context, 
since the BDL affects only collisions at central impact 
parameters, where the diffractive amplitude is anyway 
suppressed due to the vanishing rapidity gap survival 
probability (see Sec. Ivl below). 

B. Combining hard and soft interactions 

In the second step, we formalize the interplay of hard 
and soft interactions in the amplitude of the hadronic 
diffractive process, Eq. (Q. To this end, we invoke the 
parton picture of the proton wavefunction, as developed 
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soft 



FIG. 5: Schematic illustration of hard and soft interactions in the parton picture of double-gap exclusive diffractive pp scattering, 
(a) The hard scattering process producing the large-mass system (Higgs, dijet) is represented by a local operator in parton 
degrees of freedom, (b) Hard and soft interactions are approximately independent since they happen over widely different time- 
and distance scales. 



by Gribov [g. We consider the process (Q) in the CM 
frame, in which the two protons in the initial state have 
longitudinal momenta ±^/s/2, and zero transverse mo- 
mentum, Pi±,P2± = 0. Since mn <^ a/s/^i angular 
momentum conservation implies that the reaction ampli- 
tude is approximately diagonal in the transverse coor- 
dinates of the colliding protons (i.e., in impact param- 
eter) as in two-body elastic scattering at high energies. 
We thus consider partonic configurations centered around 
the transverse centers of the two protons, in which the 
partons carry fractions of the longitudinal momentum of 
the respective proton. We may regard the hard scatter- 
ing process as an operator in the basis of these partonic 
states, denoted by Vhard- Soft interactions, which build 
up the partonic wavefunctions, are governed by a soft 
Hamiltonian, i?soft- While we do not know their explicit 
form, we can nevertheless state some important proper- 
ties of these operators: 

1- Miard is local in time (instantaneous) on the typical 
timescale of soft interactions, 

2. Vjiard is local in transverse position on the distance 
scale over which the transverse spatial distribution 
of partons in the protons changes due to soft inter- 
actions, 

3. Vhard preserves the number of partons, since scat- 
tering of both gluons from the same parton domi- 
nates in the hard regime. 

4. Vhard preserves the helicity of the colliding par- 
tons, because in a perturbative gauge theory the 
dominant contribution to the interaction of partons 
over large rapidity intervals comes from the parton 
helicity-conserving component of the propagator of 
the exchanged gluon |4^ . 



Hamiltonian of soft interactions, 

[Vhard, Hsoh] = 0. 



(52) 



This is the mathematical statement of the independence 
of hard and soft interactions in diffraction. A schematic 
illustration of this picture is given in Fig. [S] We shall 
refer to Eq. (|52() and the formulas derived from it as the 
independent interaction approximation. 

The amplitude for the double-gap exclusive diffractive 
process is then determined by the matrix element 

Tdifi' = (piP2 I 'S'soft(oO,0) Vliard S'soft (0,-CX)) IP1P2), 

(53) 

where 

Ssottit2,h) = T / dt exp{itHsott) (54) 
Jti 

is the time evolution operator due to soft interactions 
(we have put the time of the hard interaction at t = 0). 
Because of Eq. the operator Vhard commutes with 

the soft time evolution operator. Using the property 

'5'soft(oO,0)S'soft(0, -Oo) = S'soft(oO, -CXj) = S'soft, 

(55) 

where S'soft is the S-matrix due to soft interactions, we 
obtain 



Tdiff = (P1P2I ^hard Ssoft \PlP2) ■ 



(56) 



Thus, the amplitude is expressed in terms of the observ- 
able matrix elements of the soft S-matrix, and the op- 
erator Vhard, calculable in QCD at the partonic level. 
Eq. (|56|l is the fundamental expression for discussing the 
physics of RGS within our approach. 



C. Suppression of inelastic diffraction 



As a consequence of these properties, we conclude that 
the operator of the hard reaction commutes with the 



In the next step, we evaluate the amplitude for 
the double-gap exclusive diffractive process based on 
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Eq. H56() . by inserting "intermediate" states (actually, 
states at t — oo) between the operators. In principle, 
one needs to sum over all diffractive states (elastic and 
inelastic) produced by the operator Vhard- An important 
question is which states can give large contributions to 
the matrix element. In fact, it turns out that the different 
preferences of hard and soft interactions severely restrict 
the range of states which can effectively contribute. 

Simple arguments show that large-mass diffractive 
states should make a negligible contribution to Eq. H56|l . 
If the two hard gluons in the hard interaction are at- 
tached to two different partons in the proton, the in- 
elastic states predominantly produced are two jets and 
gluon bremsstrahlung. It is virtually impossible to pro- 
duce such states in soft interactions, hence they cannot 
contribute to Eq. H56|l . If the two hard gluons are at- 
tached to the same parton, the cross section of inelastic 
diffraction is small for small t because of the small over- 
lap integral with the inelastic state (most of the overlap 
is with the elastic state), while for large t one produces 
a single parton with transverse momentum p± ~ \/— i, 
which again is a state difficult to reach through soft in- 
teractions. In addition, for t ^ soft diffraction at LHC 
energies is known to be dominated by the spin-flip ampli- 
tude, which further suppresses the overlap integral. To- 
gether, this restricts the possible mass range of diffrac- 
tively produced states to ~ few GeV^ . 

For a more quantitative estimate, we suppose that the 
state produced through inelastic diffraction has the form 
\pp) + ^\pX)i where the state X is different from the 
proton, and e is a small correction. We can then estimate 
e from the Schwarz inequality: 



£ _ / 0"soft(PP Xp) Crhard(PP ^ Xp) 
2 V CFsoftiPP ^ PP) CThaidiPP ^ PP) 



(57) 



Eq. H56|) into a tool for quantitative evaluation of the 
diffractive amplitude and the RGS. In particular, with 
the pp intermediate state we can approximate the matrix 
element of the soft-interaction time evolution operator by 
that of the full ^-matrix, i.e., the pp elastic scattering 
amplitude, which is known experimentally; see SectionlTTl 
For this approximation to be legitimate it is crucial that 
scattering at small impact parameters turns out to be 
strongly suppressed due to the approach to the BDL in 
pp elastic scattering, as will be seen from the results of 
Section The diffractive process is thus dominated by 
large impact parameters, where pp elastic scattering is 
dominated by soft interactions. 

In the studies of double-gap exclusive diffraction based 
on the pomeron model of Ref. inelastic intermedi- 
ate states were effectively included by way of a two- 
component formalism (however, no explicit non-diagonal 
"transition" GPDs were introduced). We have argued 
here that these contributions are strongly suppressed, be- 
cause of the small overlap of states accessible in hard and 
soft interactions. We shall comment on the implications 
of this for the numerical values of the RGS probability in 
Section EH 



D. Evaluation of the diffractive amplitude 

It remains to actually evaluate the matrix element H56(l 
with \pp) intermediate states, using the specific form of 
the hard scattering amplitude and the pp elastic scatter- 
ing amplitude. We insert a set of pp intermediate states 
in the form 



(Ppl f (fip2 



{2nY^s J (27r)3yi 



biP2>(PiP2l> (58) 



where 



where (Thard(pp ^ pX) is the cross section for the pro- 
duction of the state \pX) by the operator Vhard- Analysis 
of the Tevatron data (for a review, see Ref. 45]) shows 
that the fraction of diffractive events in soft collisions de- 
creases with increasing energy, and that the distribution 
over the excitation mass is oc 1/M|-. As a result, we 
expect that at the LHC energy e < 2 • 10"^. Thus, the 
diffractively produced state is actually the \pp) state, and 
the contributions from inelastic diffraction are small. 

The small overlap between hard and soft diffraction 
can also be understood as the result of the different im- 
pact parameter dependence of both types of processes. 
Hard diffraction occurs mostly at small impact parame- 
ters, 6^ ~ Bg. Soft diffraction, because of the approach The factor 
to the BDL, occurs mostly at large impact parameters, 
6^ ~ -B, which, moreover, rapidly grow with the colli- 
sion energy. We note again that the peripheral nature of 
soft diffraction was established already within Reggeon 
field theory, where it was found that the BDL solves the 
consistency problem associated with the triple Reggeon 
formula [20| . 

The restriction to the pp intermediate state turns 



where we have approximated the energy of the individ- 
ual protons by \fs/2. The matrix element of the opera- 
tor l^ard between the two-proton states is, by definition, 
given by [cf. Eqs. (HHJ] 

(P1P2I ^hard |PlP2> = '«(s,6,6) 

X Fg(x2,6,i2;Q'), (59) 



ii = -(Pi± -Pl_L)^ 



n{s,S.l,S.2) = Chard 

X Hg{xi,Ci,ti = 0) 

X Hg{x2,C2,t2 = 0) 



(60) 
(61) 



(62) 



represents the symbolic result for the absolute normaliza- 
tion of amplitude of the hard scattering process; it con- 
tains the amplitude of the two-gluon exchange process, 
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FIG. 6: The amplitude for double-gap exclusive hard diffraction in momentum representation, Eqs. 16611 - 1691 . The first term 
is the amplitude of the hard reaction alone, the second term the correction due to soft elastic rescattering. Only the transverse 
momenta of the protons are indicated; the momentum transfer due to soft elastic scattering is A_l. 



Chard, including the information about the ggH coupling 
given by the electroweak theory, as well as the informa- 
tion about the gluon GPD in the colliding protons at 
i = 0. The information about the transverse momentum 
dependence of the amplitude is contained in the two- 
gluon formfactors, Fg, cf. Ea. (|19|l . Furthermore, we re- 
place in Eq. (I5(: 



'5'soft 



S 



1 + (5-1), 



(63) 



and use the fact that the matrix element of the operator 
5 — 1 is given by 



{p;p;\s~i\p^P2) = *(2^)4 5(4)(p;' 

X (47r)Tei(s,t), 



'P2 -Pl 



with 



t = -(Pl±-Pl±)^ = -(P2±-P2±)^ 



P2) 

(64) 



(65) 



Finally, taking into account that at high energies the 
energy-conserving delta function in Eq. 1)64(1 effectively 
conserves longitudinal momentum, and combining the 
contributions from the two terms in Eq. H63() . we obtain 



fd^A 

7diff(Pi±,P2±) = J 7^ 



KFg {xu^iJi; g2) Fg(x2,6,^2; g') 
(27r)2<5(2)(A^) + ^re,(s,i) 



where now 



h = 

t EE 



-{P'll 
-{P'21 



A,)2, 



(66) 



(67) 
(68) 
(69) 



This result has a simple interpretation (see Fig.O. The 
first term in the brackets represents the amplitude of 
the hard reaction alone. The second term represents the 
contribution in which the hard reaction is accompanied 
by soft elastic rescattering with transverse momentum 
transfer A_l. The total amplitude is the coherent su- 
perposition of the two contributions. We note that the 
form of this result is analogous to the well-known ab- 
sorption corrections in Regge phenomenology, in which 



an elementary Regge pole amplitude is modified by elas- 
tic rescattering. 

It is instructive to convert the result (|66|l to the trans- 
verse coordinate representation. Substituting the Fourier 
representation of the gluon GPDs, Eq. (|22|l . and the 
representation of the pp elastic amplitude in terms of 
the profile function, Eq. and using standard Fourier 
transform manipulations, we obtain 



d b dpi dp2 



'7diff(p'i_L,P2±) ^ 

X K Fg {xi,£,i,pi; Q'^) Fg {x2,^2,P2; Q'^) 

X [i-r(s,6)]. 



(70) 



Here the scattering amplitude is represented as the co- 
herent superposition of amplitudes corresponding to pp 
scattering at given transverse displacement (impact pa- 
rameter), b. The amplitude for the hard process is pro- 
portional to the product of the transverse spatial gluon 
transition densities at positions pi^2 relative to the cen- 
ters of the respective protons, with the three transverse 
vectors satisfying the triangular condition pi — p2 — b 
(see Fig.[71l. The modifications due to elastic rescattering 
now take the form of a multiplication of the hard scatter- 
ing amplitude with the "absorption factor," 1 — T{s,b). 
Note that the modulus squared of this factor can be in- 
terpreted as the probability for "no inelastic interaction" 
in pp scattering at a given impact parameter, cf. Eq. 10}. 
This interpretation will be explored further in Sec. 

Our partonic approach allows us to calculate the am- 
plitude for double-gap exclusive diffraction in a model- 
independent way in terms of the gluon GPD and the 
phenomcnological pp elastic scattering amplitude; see 
Eqs. jnnHZni- R-ef. SJ such processes were studied 
using a model of elastic pp scattering which included the 
enhanced eikonal series of single Pomeron exchanges and 
the triple-Pomeron vertex to describe the soft specta- 
tor interactions. The expression for the amplitude in the 
case of a single Pomeron exchange in that model (and 
without inelastic intermediate states) would formally co- 
incide with our expressions (|66H70I) . Whether the same 
is true for the full amplitude in that model is less clear; 
cf. the discussion of the numerical results in Section [VBI 
below. 
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Hard 
process 




FIG. 7: Illustration of the transverse coordinate representa- 
tion of the diffractive amplitude, Eq. 17011 . The hard scatter- 
ing process is local in transverse space. The centers of the 
colliding protons are displaced by the distance b — \b\, and 
Pi. 2 ~ IpiM are the distances from the centers to the point 
of the hard process. 



THE RAPIDITY GAP SURVIVAL 
PROBABILITY 



We now use our general result for the amplitude of 
double-gap exclusive diffractive processes in the inde- 
pendent interaction approximation, Eq. H7U|) . to calculate 
the RGS probability for such processes. At this level of 
approximation, we shall recover a simple "geometric" ex- 
pression for the RGS probability, which permits a prob- 
abilistic interpretation and was heuristically derived in 
Refs. 0, 0, 0. We discuss the impact parameter de- 
pendence of RGS and stress the the crucial role of the 
BDL in stabilizing the numerical estimates ensuring a 
model-independent result. 



A. Impact parameter representation 

In order to compute the cross section for double-gap 
exclusive diffractive production of a given state at fixed 
rapidity, we integrate the modulus squared of the ampli- 
tude IZOI) 

over the final proton transverse momenta. By 
standard Fourier transform manipulations we obtain 



cTdiff = (kinematic factors) x / - — / 

J V'T^r J 



(2^)2 



X \TMp'i^,P2i^)f (71) 
= (const) X J d?p^ J d?p2 F^{p,) F^ip2) 



X |i-r(p2-Pi)r 



(72) 



(for brevity we suppress all arguments except the trans- 
verse coordinates in Fg and F). The RGS probability due 
to soft interactions 0, by definition, is given by the ratio 
of the cross section of the physical double-gap diffractive 
process (|72|) to the cross section of the hypothetical pro- 
cess with the same hard scattering subprocess but with 
no soft spectator interactions, corresponding to expres- 



sion lO with F = 0, 

^2 _ o-diff (physical) 



f diff (no soft interactions) 



(73) 



We can cast this ratio in a simple form. We rewrite the 
convolution integral in Eq. H72I) by inserting unity in the 
form (c/. Fig. CI) 



J dHS^^\b + pi- p2), 



(74) 



and introduce a normalized impact parameter distribu- 
tion, 

Phard(b) = Jd^Pl j P2 S^^Hb + - P2) 



F'iP2) 



J d'p[ F^{p[) J d'p'2 F^{p'2) 



(75) 



which satisfies 



J dHPi,,,d{b) = 1. 



(76) 



In terms of this distribution the RGS probability l|73l) is 
expressed as 



d'foPhard(b) |l-F(b)|^ 



(77) 



This result agrees with the expression for the RGS prob- 
ability derived heuristically in Refs. 0, ^ • For 
the comparison of our result for the RGS probability with 
that obtained with the pomeron model of Ref. Q we re- 
fer to Section FV Bl below: see also the comments at the 
end of Sec. HVDl 

Expression (|77|l for the RGS probability allows for a 
simple probabilistic interpretation. Consider a pp colli- 
sion at given impact parameter, b — \b\. Since the hard 
two-gluon exchange process is effectively local in trans- 
verse space, the probability for it to happen is propor- 
tional to the product of the squared transverse spatial 
distributions of gluons in the two colliding protons, in- 
tegrated over the transverse plane, as given by the nu- 
merator of Eq. H75|l . Consider now a hypothetical sample 
of pp events with the two-gluon induced hard scattering 
process, but an otherwise arbitrary (non-diffractive) fi- 
nal state. By the laws of probability, the distribution of 
impact parameters in this sample is given by the normal- 
ized distribution Phard(b), Eq- lf75|l . A diffractive event 
results if the spectator systems of the two protons do not 
interact inelastically. The probability for this to happen 
in a pp collision at fixed b is given by |1 — F(b)p, cf. 
Eq. , in analogy to the well-known formula for inelas- 
tic scattering in non-relativistic theory |4^. The RGS 
probability, which is defined as the fraction of diffractive 
events in the sample of all events containing the same 
hard scattering process, is then given by the average of 
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FIG. 8: The integrand (impact parameter distribution) in 
the RGS probability, Eq. H77|l . for Higgs boson produc- 
tion at the LHC energy. Dashed hne: 6-distribution of 
the hard two-gluon exchange, Phard(6), Eq. 17511 . evaluated 
with the exponential parametrization of the two-gluon form- 
factor, Eq. (Oni with Bg = 3.24 GeV"^. Solid line: The 
product Phard(&)|l — T{b)\^, evaluated with the exponential 
parametrization, Eq. 11211 . with B — 21.8 GeV~^. The van- 
ishing of |1 — r(b)p, at small b, cf. Fig.Q strongly suppresses 
contributions from small impact parameters. Note that the 
plot shows 2nb times the functions of impact parameter; the 
small-& part of the dashed curve [distribution Phard(&)] would 
be close to the left boundary of the plot and was omitted for 
better legibility. The RGS probability, S^, is given by the 
area under the solid curve. 



vanishes for <C 1/B where pp scattering approaches 
the BDL. The product of the two probabilities is sup- 
pressed both at small and at large b, and thus concen- 
trated at intermediate values of b. 

This point can nicely be illustrated with the Gaussian 
parametrizations of the transverse spatial distribution 
of gluons, Eq. (|34|l . and the pp elastic profile function, 
Eq. (|12|) . With the Gaussian form (|34|l . the convolution 
integral in Eq. (|75|l can be computed analytically, 

_ exp[-bV(2i3,)] 
^''^''^^^^ - 2^, ■ 

This function is shown by the dashed line in Fig. |H1 The 
integrand of Eq. H77|) is given by 

Phard(b) |l-r(b)|^ 

and is shown by the solid line in Fig. |S1 It is suppressed 
both for 6^ <C 1/i? (because of the "blackness" of the pp 
amplitude) and for &^ 1 /Bg (because of the vanishing 
of the overlap of the two gluon distributions), and thus 
concentrated at intermediate values of b. The maximum 
of 27r& times the combined distribution is at 

b^ w 5Bg {Bg<^B). (80) 



this function with the normalized 6-distribution in the 
sample, Eq. lf77|l . 

It needs to be stressed that the impact parameter of 
a single pp event is not observable, being a microscopic 
quantity beyond the reach of any experimental appara- 
tus. In the above arguments, the impact parameter plays 
the role of a randomly chosen external parameter. How- 
ever, using information about the transverse spatial dis- 
tribution of gluons in the proton from independent mea- 
surements, we can calculate the probability for certain 
hard processes in a pp collision as a function of the im- 
pact parameter, and thus infer the distribution of impact 
parameters in a sample of events with the same hard pro- 
cess. This logic was used in Ref. to devise a trigger 
on central collisions in inclusive pp scattering by requir- 
ing hard dijet production at small rapidities. Here we use 
the same strategy to model soft spectator interactions in 
double-gap exclusive diffr active pp scattering. 

The integrand in Eq. H77(l describes the effective distri- 
bution of impact parameters in a sample of double-gap 
diffractive events, and reflects the interplay of hard and 
soft interactions at the cross section level. The proba- 
bility for the hard process, -Phard(b), favors small impact 
parameters, which maximize the overlap of the large- 
X gluon distributions in the protons, and vanishes for 
b^ ^ 1/Bg. The probability for no inelastic soft interac- 
tions, |1 — r(6)p, favors large impact parameters, which 
increase the chances for the protons to stay intact, and 



We see that within our two-scale picture of the trans- 
verse structure of hard and soft interactions, cf. Fig. |21 
the dominant impact parameters in double-gap exclusive 
diffractive processes are determined by Bg — the smaller 
of the two areas — , but may be numerically large be- 
cause of a large numerical factor. The RGS probability, 
Eq. (|77|l . is given by the integral of 27r& times Eq. (|79() 
{i.e., the area under the solid curve in Fig. ISJ, and can 
be computed analytically, 

2 2B^g 
^ ^ {B + Bg){B + 2Bg) -W (^^«^)- 

(81) 

The gap survival probability is of the order {Bg/B)^, 
i.e., it is proportional to the square of the ratio of the 
transverse area occupied by hard gluons to the area cor- 
responding to soft interactions. Thus, our two-scale pic- 
ture offers a parametric argument for the smallness of the 
rapidity gap survival probability. 

The approach to the black-disk limit in pp scattering 
at high energies, i.e., the fact that r(6) 1 at small b, 
plays a crucial role in determining the numerical value of 
the RGS probability and ensuring stability of our calcula- 
tion with respect to variation of the parameters. A small 
deviation of the profile function from unity at 6 = 0, of 
the form r(b = 0) = 1 — e with e <C 1, would change the 
result for the gap survival probability to 

^ S\uL + (82) 
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FIG. 9: The RGS probability for double-gap exclusive diffrac- 
tive processes, Eq. 1771 . as a function of the squared CM 
energy, s. The Tevatron and LHC energies are marked by 
arrows. Shown are the results obtained with the dipole 
parametrization of the two-gluon formfactor 13111 , for different 
values of the mass parameter, m'g. The value rrig — 1 GeV'^ is 
appropriate for Higgs boson production at the LHC at central 
rapidities. The profile function of the pp elastic amplitude was 
taken from Ref. [l^ (c/. Fig.[TJ. 



FIG. 10: Dependence of the RGS probability for double-gap 
exclusive diffractive Higgs production at the LHC, Eq. 177^ . 
on the mass parameter in the dipole parametrization of 
the two-gluon formfactor, m^, Eq. 13111 . Also shown are 
the results obtained with the exponential parameterization, 
Eq. lESJ, with Bg = 3.24/m^, cf. Eq. ^ and Fig. El The 
profile function of the pp elastic amplitude in both cases is 
the one of Ref. |l3|. 



[here we have used that Bg ^ B, and that the integral 
of -Phard is Unity, cf. Eq. (|7^ ]. The approach to the BDL 
effectively eliminates r(6 = 0) as a parameter in our cal- 
culation. We stress again that the experimental data as 
well as theoretical arguments indicate that the BDL is in- 
deed reached in pp scattering at small impact parameters 
at the LHC energy. 



B. Numerical estimates 

For a numerical estimate of the gap survival probability 
we evaluate Eq. H77|l with the dipole parametrization of 
the two-gluon formfactor, Eq. H31|l . and the parametriza- 
tion of the pp elastic amplitude of Ref. For Higgs 

production at the LHC {^/s = 14TeV) at central rapidi- 
ties the momentum fractions of the annihilating giuons 
are xi^2 ~ 10~^ (at a scale <C m'jj). For such val- 
ues of X in principle the contributions of the nucleon's 
pion cloud to the gluon density at transverse distances 
p ~ 1/{2M^) need to be taken into account; see Sec- 
tion ^O] As we shall explain below, these contributions 
to the gluon density involve correlations in the nuclcon 
wavefunction, which effectively reduce their contribution 
to RGS, and should not be included in the estimate based 
on Eq. (|77|) . We therefore use in our estimate at the LHC 
energy the simple dipole formfactor with w 1 GeV^ , 
which does not include the pion cloud contribution. With 
this choice of parameters Eq. I77|) gives for the RGS prob- 
ability for Higgs production at the LHC 



The energy dependence of the RGS probability is shown 
in Fig. O for various values of the mass parameter of 
the two-gluon formfactor, m^. At the Tevatron energy 
{y/s = 1.9 TeV), the gluon momentum fractions x\;i are 
of the order ^ 10^^, for which the pion cloud contribu- 
tions to the gluon density are naturally absent. While a 
mass parameter = 1 GeV^ is still reasonable in this 
situation, even higher values might be relevant in this 
case. Taking into account this effective change of 
with s via the xi, a;2-dependence, the actual variation 
of the RGS probability between the LHC and the Teva- 
tron energies is less pronounced than it appears from the 
fixed-TTig curves in Fig. |51 

FigurelTUlfsolid line) shows the dependence of the RGS 
probability at the LHC energy on the mass parameter in 
the dipole parametrization of the two-gluon formfactor, 
m^, Eq. The value of S*^ strongly increases with 

decreasing m^, i.e., with increasing radius of the the 
transverse spatial distribution of giuons, similar to the 
behavior found with the simple exponential parametriza- 
tions, Eq. (|81|l . To illustrate the sensitivity of our numer- 
ical predictions to the shape of the two-gluon formfactor, 
we also show the results for S"^ obtained when replacing 
the dipole formfactor by the exponential, Eq. 1)33(1 . with 
Bg = 3.24:/ml, cf. Eq. (Figure [TOI dashed line). 

One sees that the numerical values are rather different, 
in spite of the apparent similarity of the two formfactors 
over a wide range, \t\ < 1 GeV^ {cf. Fig.EJ. Comparing 
the two curves of Figure [TUI we estimate the uncertainty 
of our numerical prediction for the RGS probability due 
to the uncertainty of the shape of the two-gluon form- 
factor to be at least ~ 30%. 



= 0.027. 



(83) 



In a similar way, we can estimate the uncertainty of 
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the RGS probability due to the uncertainty of the profile 
of the pp elastic amplitude. Comparing the numerical 
values obtained from Eq. (|77|l with the parametrization 
of Ref. Jit) and with the simple exponential parametriza- 
tion, Eqs. Hll|l and H12() . for the same two-gluon form- 
factor, we again find variations of the order of ~ 30%. 

The relatively high sensitivity of the numerical esti- 
mates to the shape of the two-gluon formfactor and the 
profile function of the pp elastic amplitude can be under- 
stood as a result of the peculiar interplay of hard and soft 
interactions in Eq. (|77() . The different impact parameter 
dependence of hard and soft interactions (c/. Fig. IHland 
the discussion above) essentially eliminates contributions 
from the regions corresponding to the "bulk" of the in- 
dividual distributions, Pha,id{b) and |1 — r(6)|, allowing 
for significant strength only in an intermediate region of 
impact parameters, where there is considerable sensitiv- 
ity to the shape of the distributions. This seems to be 
a principal feature of estimates of the RGS probability 
based on Eq. ifTTji . 

Detailed numerical studies of the RGS probability were 
made within the eikonalized Pomeron model for soft in- 
teractions, see Ref. for a summary. We would like 
to comment on these estimates from the perspective of 
our approach. As we already noted, the approach of 
Ref. |j| includes contributions from inelastic intermediate 
states (albeit without introducing non-diagonal "transi- 
tion" GPDs). We have argued that within the indepen- 
dent interaction approximation these contributions are 
strongly suppressed and should not be included, see Sec- 
tion IIV CI Nevertheless, the numerical result for the 
RGS probability in Higgs boson production in double- 
gap diflfraction at the LHC quoted in Ref. ^, ^ 0.023 
|66l |. is rather similar to our estimate H83|l . Note that 
Ref. Q ascribes an uncertainty of ~ 50% to this value; 
we have estimated a similar uncertainty for our result 
due to the combined uncertainties in the profile function 
and the two-gluon formfactor (see above). It is interest- 
ing to ask why the results are so similar when the two 
approaches differ in their treatment of inelastic diffrac- 
tion. In order to clarify this question, we have evalu- 
ated our expression for the RGS probability, Eq. (f77|) . 
with the profile function of the pp elastic amplitude ob- 
tained within the model of Ref. which was kindly 
provided to us by M. Ryskin; we emphasize that this 
is not the same as evaluating the expression for the 
RGS probability given in Ref. Using the exponen- 
tial form of the two-gluon formfactor H33(l with param- 
eters Bg = (4, 5.5, 10.1) GeV~^ we obtain in this way 
52 = (0.042, 0.069, 0.157), which should be compared 
to the results quoted in Ref. :4], = (0.02, 0.04, 0.11). 
One sees that our Eq. H77|l gives systematically larger 
values than the approach of Ref. |j| for the same pro- 
file function and the same two-gluon formfactor. This 
difference should be attributed to the effect of inelastic 
diffraction. What is then quoted as the final estimate 
of depends on the preferred value of Bg. The RGS 
probability strongly decreases with Bg, cf. Eq. H81|l and 



Fig. ^] For Higgs production at the LHC we use a value 
of Bg = 3.24 GeV~^ (corresponding to a mass parameter 
in the dipole parametrization of = 1 GeV^), which is 
based on analysis of the J/tp photoproduction data over 
a wide range of energies and takes into account the effects 
of QCD evolution {cf. Section UTTjl . With this value we 
obtain S'^ = 0.030 for the exponential two-gluon form- 
factor and the profile function of Ref. Jj. This value of 
Bg is lower than the range of values considered in Ref. Q 
(our Bg = 2b in the notation of that paper). The value 
of Bg for Higgs production at the LHC taken in Ref. ^ 
is Bg = 4GeV"^, which resuhs in 5*^ 0.02 in their 
model. One sees that the different values of Bg partly 
compensate the differences due to the treatment of in- 
elastic diffraction in the two approaches. We thus con- 
clude that the similarity of the final numerical estimate 
of Ref. (3 with our results is somewhat accidental. In 
any case, the differences between the numerical results of 
both approaches are within the estimated uncertainties. 

Potentially more important than the uncertainties of 
our calculation of the RGS probability within the inde- 
pendent interaction approximation are effects of possible 
correlations between hard and soft interactions. These 
effects can naturally be incorporated into our partonic 
picture, and further decrease the RGS probability com- 
pared to the independent interaction approximation (see 
Section rvTBl. 



VI. BEYOND THE INDEPENDENT 
INTERACTION APPROXIMATION 

Our treatment of RGS so far was based on the idea that 
hard and soft interactions are approximately independent 
because they proceed over widely different time- and dis- 
tant scales. It is clear that this approximation has certain 
limitations, concerning both the range of its applicability 
and its accuracy. In this Section we discuss various phys- 
ical effects which violate the assumption of independence 
of hard and soft interactions and give rise to corrections 
to the estimates of the RGS probability of Section 
These are (a) the increase of hard screening corrections 
with energy, (b) correlations between hard and soft in- 
teractions. These corrections have not been considered 
in previous treatments of RGS in Refs. 0, HE Ell- 

A. Hard screening corrections 

The independent interaction approximation relies on 
the assumption of widely different characteristic scales of 
hard and soft processes. However, the difference between 
the two scales tends to decrease, and may even disappear, 
at high energies, because of the fast increase of the am- 
plitudes of hard processes with energy. One example of 
this effect is the BDL in pp elastic scattering at central 
impact parameters, which can be explained both on the 
basis of hard and soft interactions (see Section • In 
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needs to follow the space-time evolution of the produc- 
tion of small-x, large-virtuality partons, and must not 
restrict the discussion to the gluon GPDs in the individ- 
ual protons. 

We emphasize that the screening effect described here 
is not included in the definition of the RGS probability 
due to soft interactions, as it corresponds to a modifi- 
cation of the probability of finding two gluons at small 
transverse distances without reference to the spectator 
interactions. Detailed calculations of the RGS probabil- 
ity including this eff'ect require modeling of the color con- 
nections in the parton wavcfunctions of the protons. An- 
other perturbative screening correction to the hard pro- 
cess, due to diffraction into high-mass states (with the 
mass increasing with the collision energy) was considered 
in Ref. (49j , which leads to a complementary suppression 
of the amplitude for exclusive double-gap diffraction. 



FIG. 11: A typical diagram describing hard screening correc- 
tions to the hard amplitude for double-gap diffractive pro- 
duction. Shown is only the generic structure of the partonic 
ladders; the dominant contribution comes from gluons. 

diffractive scattering, the increase of hard amplitudes at 
high energies leads to a reduction of the RGS probabihty 
relative to the estimates presented in Section Ivl 

One specific mechanism which can lead to qualitative 
changes of the hard scattering process at high energies is 
"local" absorption of the hard scattering amplitude for 
diffractive production. This corresponds to the familiar 
attenuation of the hard gluon exchange diagram for Higgs 
production, but by a two-gluon ladder (see Fig. lll|) . Here 
the two-gluon ladder is attached to the partons which 
emitted the gluons involved in the Higgs production. 
Actually, in order to regularize the infrared divergences 
present in the two-gluon ladder we may consider instead 
the amplitude for the scattering of two colorless dipoles, 
in which the typical virtuality of the "constituents" is k'^ ; 
a second parton with such virtuality is anyway present 
as a result of QCD evolution. The amplitude of the addi- 
tional interaction is then suppressed relative to the origi- 
nal hard amplitude by a factor al{kj_/k'jJ{xo/x)^, where 
k± is the characteristic gluon transverse momentum in 
the absorptive ladder. The characteristic fc^ increases 
with energy, and A « 0.2 - 0.25 at = 4 Q^y^ 
Ultimately, the amplitude of the additional interaction 
would thus reach a strength comparable to the maximal 
one (BDL), resulting in complete suppression of RGS. 

A semi-realistic estimate, using the eikonal approxi- 
mation for the scattering of the colorless dipoles and as- 
suming A « 0.2 — 0.25, suggests that at the LHC energy 
the additional suppression is significant (a factor 2-3) for 
Higgs boson production by valence gluons, but much less 
so for processes initiated by valence quarks; at the Teva- 
tron energy the effect should be much less pronounced 
for both gluons and quarks. In order to evaluate this ad- 
ditional suppression at LHC energies quantitatively, one 



B. Correlations between hard and soft interactions 

In our treatment of RGS in Section IIVBI we have not 
taken into account effects due to correlations between 
partons in the wavcfunctions of the colliding protons. 
Generally speaking, in the presence of such correlations 
the selection of particular configurations by the hard 
scattering process changes soft interactions as compared 
to "average" configurations. The neglect of this change is 
basic to the approximation H52I) . in which hard and soft 
interactions are assumed to commute also with respect 
to transverse degrees of freedom. A consistent treatment 
of correlation effects requires a much more detailed de- 
scription of hard and soft interactions and the parton- 
hadron interface than the one given in Section Hvl Here 
we would like to discuss this problem at the qualitative 
level, preparing the ground for a future in-depth investi- 
gation. 

As a pedagogical example illustrating the effect of cor- 
relations on the RGS probability, let us consider double- 
gap diffractive proton-deuteron (pd) scattering — pos- 
sibly including the quasi-elastic channel — within the 
framework of Eq. H77|l . The profile function for pd elastic 
scattering is significantly smaller than that for pp scat- 
tering even for energies at which pp scattering is close to 
the black-disk limit, 

Y^p ^ < 1. (84) 

This follows from Eq. when noting that the dp total 
cross section is approximately equal to twice the pp cross 
section, while the spatial size of the deuteron is much 
larger than the proton radius, and can be understood 
simply as the result of a "dilution" of the blackness of the 
individual protons due to their transverse motion in the 
deuteron bound state. At the same time, the transverse 
spatial distribution of gluons is now characterized by the 
transverse size of the deuteron. Because of Eq. (|84|) . the 
factor |1 — F'^Pp in Eq. H77|) is always of order unity, 
and no significant suppression takes place as it does in 
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pp scattering. Thus, one would conclude that the RGS 
probability is much larger in dp than in pp scattering, 
which is clearly nonsensical. The paradox is resolved by 
noting that hard and soft interactions in diffractive pd 
scattering are highly correlated. By considering events 
with a hard process one is effectively selecting configu- 
rations in which the projectile proton scatters from one 
of the nucleons. Soft interactions in these configurations 
are significantly larger than in pd scattering in average 
configurations, in which there is a substantial chance of 
the projectile "missing" the nucleons in the deuteron [6^ . 

The deuteron example shows that transverse correla- 
tions in the wavefunctions can qualitatively change the 
picture of RGS. In particular, positive spatial correlations 
between hard partons and the opacity for soft interac- 
tions decrease the RGS probability compared to the un- 
correlated estimate based on Eq. l(77|) . In connection with 
pp scattering, Eq. QJ, we now discuss the effects of two 
types of transverse correlations: (a) long-distance cor- 
relations due to scattering from the proton's long-range 
pion field ("pion cloud"), (b) short -distance correlations 
related to parton clustering in "constituent quarks" . 

A distinctive contribution to diffractive pp scattering 
at small results from the process in which a soft 
pion, emitted and absorbed by proton 1, scatters diffrac- 
tively from proton 2. This contribution could properly be 
calculated using the known pion-nuclcon coupling, and 
applying Eqs. H66() - (|69|l to the np diffractive amplitude, 
with the two-gluon formfactor in the pion and the pro- 
file function of np elastic scattering. It is generally small, 
for two reasons. First, the coupling of the soft pion to 
the nucleon is small because it is the Goldstone boson 
of spontaneously broken chiral symmetry. Second, soft- 
ness of the pion implies that its longitudinal momentum 
fraction in the proton be small, y m-j^/mp. In Higgs 
production at the LHC, where xi^2 ^ lO^'^, this puts 
the momentum fraction of the gluons in the pion at rel- 
atively large values, z ~ xi^2/y ^ 10^^, where the gluon 
distribution is not enhanced by DGLAP evolution. 

In the partonic picture, the pion cloud contribution 
represents the result of specific correlations between hard 
and soft partons in the proton wavefunction. This obser- 
vation has an interesting consequence for the estimate of 
the RGS probability for pp diffractive scattering based 
on Eq. H77|l . Namely, the gluon distribution in the pro- 
ton receives a contribution from the pion cloud at trans- 
verse distances p 1/(2/71^) and momentum fractions 
X ^ niTj/mN [33 ■ Including this contribution in Eq. I|77l) 
would be inconsistent, since a hard process involving 
these gluons in the pion cloud should be accompanied by 
a very specific modification of the soft interactions, which 
is not accounted for in Eq. (|77|l . For this reason, we have 
not included an explicit pion cloud contribution in the 
two-gluon formfactor parametrization used in our esti- 
mate leading to the value (|83l) . More precisely, if we knew 
the "physical" gluon GPD, which by definition includes 
the pion cloud contribution, we would need to remove 
this contribution before using the GPD in Eq. (|77f) . and 
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FIG. 12: The two-component parametrization of the trans- 
verse spatial distribution of gluons including the pion cloud, 
Eqs. and JHSJ. Shown are the radial distributions 

Fg{p), the integral of which determines the average glu- 
onic transverse size, (p^). For the sake of comparison, the 
figure shows the combined (bulk -I- cloud) distribution with- 
out adjustment of the normalization factor (A'' = 1), so that 
the bulk contribution is the same as in the case of no pion 
cloud. 



thus obtain a lower value for the RGS probability. This is 
just another example of the general rule that correlations 
lower the RGS probability compared to the independent 
interaction approximation, Eq. H77|l . 

To estimate the correction resulting from the removal 
of the pion cloud contribution from the gluon GPD, we 
evaluate Eq. (|77|l with a "two-component" parametriza- 
tion of the two-gluon formfactor {viz. the transverse spa- 
tial distribution of gluons) including the pion cloud. We 
write 



Fg{p) = A^[Fbuik(p) + Fcioud(p)] 



(85) 



where A'^ is a factor which ensures overall normalization 
to / d'^pFg{p) = 1. The "bulk" contribution to the two- 
gluon formfactor we parametrize by the dipole formfactor 
Eq. H31I) . the "cloud" contribution as 



with TO^ = 1 GeV^ 



^cloud (p) = 



loud 



p2 g-2m,p 



+ Pq 2m-^p 



(86) 



This form is essentially the asymptotic form of the gluon 
density at p ^ l/(2m^) for x ^ m^/mN ( "Yukawa tail" ) 
[3^ . regularized at small p such as to avoid a large con- 
tribution in the bulk region; the parameter pp is chosen 
of the order (p^)g,buik = = 0.3 fm^. The coefficient 

C'cioud we determine such that the inclusion of the cloud 
contribution increases the overall gluonic transverse size 
of the proton by 30% , which is the value found in the 
calculation of Ref. [s^l (based on the phenomenological 
gluon distribution in the pion) and supported by the J ftp 
photoproduction data; see Section UTTl Figure IT^ shows 
the transverse spatial distributions for "bulk only" and 
"bulk + cloud" , multiplied by 27rp'^, the integral of which 
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FIG. 13: The integrand (impact parameter distribution) in 
the rapidity gap survival probability, Eq. 1)77^ . for the two- 
component form of the transverse spatial distribution of glu- 
ons. The parameters are the same as in Fig. 1121 



determines {p'^)g. Figure IT!^ shows the impact parameter 
distributions in the RGS probabihty (c/. Fig. IS)) obtained 
in the two cases. One sees that removal of the pion cloud 
contribution indeed reduces the RGS probability. The 
numerical effect turns out to be rather small, 



5*2 (bulk only) 
52 (bulk + cloud) 



0.94. 



(87) 



This can be explained by the fact that the pion cloud 
contribution to the gluon density is noticeable compared 
to the bulk only at transverse distances p ^ 1/(2to^) — 
0.7 fm, while the RGS integral is dominated by rather 
short distances, peff ~ ~ 0.4 fm, c/ Figs. QandlHl 

Corrections to the rapidity gap survival probability as 
given by Eq. H77|l can also come from transverse short- 
distance correlations in the proton wavefunction (corre- 
lation length <C proton size). The Tevatron CDF data on 
pp collisions with multiple hard processes (three jet plus 
photon production) 12] indicate the presence of signifi- 
cant correlations between the transverse positions of hard 
partons over distances r ~ 0.3 fm Given one hard 

parton with x > 0.05 at a certain transverse position, it 
is much more likely to find in the proton wavefunction 
a second hard parton x ^ 0.05 within a distance ~ r 
than elsewhere in the transverse plane. As a result, in 
pp events with (at least) one hard process the probability 
for a second hard interaction is substantially larger (by 
a factor of ^ 2) than it would be without correlations. 
One may also suppose that the local cross section density 
(opacity) for soft inelastic interactions is higher near the 
position of a hard parton than elsewhere. Such corre- 
lations would result in a higher probability for inelastic 
interactions in pp events with a hard process (such as the 
hard two-gluon exchange in double-gap diffractive pro- 
duction) as compared to generic pp collisions, and would 
thus decrease the RGS probability compared to Eq. (|77|l . 

The corrections to Eq. (|77|l due to short-distance cor- 
relations can be viewed as an effective reduction of the 



size of the diffractively scattering system, from the pro- 
ton radius to the size of the transverse correlation, r. 
The corrections could be particularly large in the situa- 
tion when the correlated areas are "black spots," while 
the proton overall is still "gray" because of the dilution 
by the transverse motion. This situation would in a sense 
correspond to the above example of the deuteron, with 
the proton now playing the role of the deuteron. A quan- 
titative description of these effects would require detailed 
modeling of the correlation between hard partons and the 
opacity for soft inelastic interactions, including an anal- 
ysis of pp elastic scattering allowing local fluctuations in 
opacity, which are outside of the scope of the present 
paper. 

A model of the proton accounting for short-distance 
correlations between partons is the so-called chiral 
quark-soliton model |5^, which describes the proton 
as a system of constituent quarks bound by a classical 
pion field, see Ref. [53i] for a review. This model imple- 
ments the short-distance scale related to the spontaneous 
breaking of chiral symmetry, which appears here as the 
"size" of the constituent quark, r ^ 0.3 fm (in the Eu- 
clidean formulation of QCD this scale can be associated 
with the average instanton size in the vacuum). This 
model provides a consistent description of the twist-2 
quark- and antiquark distributions in the nucleon at the 



scale p 



600 MeV |5^]. It also suggests that the 



gluons at this scale are "packaged" inside the constituent 
quarks and antiquarks 53]. Assuming perturbative QCD 
evolution to be applicable starting from this scale one 
would thus expect significant correlations between the 
positions of hard quarks and gluons and the opacity for 
soft interactions over a transverse size r ~ 0.3 fm. In- 
corporating these correlations into the description of the 
BDL in high-energy pp scattering and the theory of RGS 
in diffractive processes is an important problem for future 
studies. 

To summarize, the dynamical mechanisms which we 
have discussed here (and indeed all mechanism which we 
are aware of) give rise to positive correlations between 
the transverse position of hard partons and the opacity 
for inelastic interactions. We can thus say with some 
confidence that estimates based on the independent in- 
teraction approximation, Eq. H77|l . represent an upper 
bound on the RGS probability. 



VII. TRANSVERSE MOMENTUM 
DEPENDENCE 

The interplay of hard and soft interactions in exclusive 
double-gap diffraction not only causes the suppression of 
the integrated cross section summarized in the RGS prob- 
ability, but also gives rise to a distinctive dependence of 
the cross section on the final proton transverse momenta. 
By observing this "diffraction pattern" one can perform 
detailed tests of the diffractive reaction mechanism, and 
even extract information about the two-gluon formfac- 
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tors of the colliding protons. The transverse momentum 
dependence also contains information about the quan- 
tum numbers (parity) of the produced system We 
consider here production of a 0"*" system, for which the 
hard scattering amplitude depends on the proton trans- 
verse momenta only through \p'ij_\ and |p2j^|, and the 
diffractive amplitude is given by Eq. (j^^; we comment 
on production of 0~ systems at the end of this section. 

For a quick orientation over the transverse momen- 
tum dependence, we evaluate the amplitude with the 
exponential parametrizations of the two-gluon formfac- 
tor, Eqs. H33|l . and the pp elastic scattering amplitude, 
Eq. (|12|) . In this case the convolution integral in Eq. Hfififl 
reduces to a Gaussian integral, and we obtain a closed 
expression for the amplitude, 



7diff(pi_L,P2J 



K exp 



X 1 



B 



Bt. 



■ exp 



2 2 1 



2Bt. 



where Bgi = Bg(xi) and Bg2 = Bg{x2) are the slopes 



corresponding to the momentum fractions xi,2 
and 



Bt, 



Bql + B, 



S2 



B. 



6:2, 



(89) 



[As a check, we note that the integral of the square of this 
expression, divided by the corresponding expression for 
5 = (no soft interactions) , reproduces the result for the 
RGS probability obtained in the coordinate space calcu- 
lation, Eq. (|HT|l .] The amplitude (|5S|l vanishes trivially 
for large \pii] or |p2_l|, independently of the directions 
of the momentum vectors. In addition, it has a zero at 
finite values of the transverse momenta, namely at 

iBglp\^-Bg2p'2,_f 

Btot 



2Btot In 
2(^91 



B 
Ba2) 



(B»Bgl,Bg2). (90) 



This zero arises because of the destructive interference 
of the amplitude of the hard scattering process alone 
and the amplitude including soft elastic rescattering, cf. 
Fig. |S1 and directly reflects the interplay of hard and soft 
interactions. It leads to a dip in the diffractive cross 
section, and thus to a typical "diffraction pattern" in 
the dependence on the transverse momentum of the first 
proton, p'l^, at fixed transverse momentum of the sec- 
ond proton, P2±- Figure H"!] shows this diffraction pat- 
tern in the kinematics of Higgs production at the LHC 
at zero rapidity, for which Bgi — Bg2. One sees that 
the pattern in P2j_ evolves from a rotationally symmetric 
one for p\j^ = to a two-centered one for large 
This basic feature of the diffractive cross section does 
not depend on the details of the parametrization of the 



two-gluon formfactor; similar results are obtained with 
the dipole parametrization, cf. the detailed comparison 
below. 

The diffraction pattern of Fig. E| implies a strong 
angular dependence of the cross section, which, more- 
over, changes with the magnitude of the transverse mo- 
menta, \p'ij_\ and \P2±\- This is illustrated in Figure [T^ 
which shows the dependence of |Tdiffp on the angle be- 
tween the transverse momenta, 0, for various values of 



IP 



\P'2J 



For small values of the transverse mo- 



menta the cross section is maximal at zero angle; for 
large values (where the cross sections as a function of 
angle runs through the diffractive dip) it is maximal at 
(f) = TT. This dependence needs to be taken into account 
when attempting to maximize the diffractive cross sec- 
tion in the search for new particles. 

An interesting questions is which specific features of 
the transverse momentum dependence of the diffractive 
cross section could be used to test the diffractive reac- 
tion mechanism, and possibly extract information about 
the two-gluon formfactors of the colliding protons. Such 
studies would be feasible in diffractive dijet production, 
which has a relatively large cross section and also allows 
one to vary the invariant mass of the diffractively pro- 
duced system, and thus the effective values of the gluon 
momentum fraction, xi and X2- To address this question, 
we again start from the explicit expression for the am- 
plitude obtained with the exponential parametrization 
of the two-gluon formfactor, Eq. H88|) . We rewrite it in 
terms of the center-of-mass and relative transverse mo- 
menta of the final-state protons, 



fP2±)/2, 



(91) 



and obtain 



7diff(p'i±,P2±) 



K exp 



Bgi + Bg2 



B 



Bu 



■ exp 



{Bgi ~ Bg2)^ p2 
2Btr,t 



{Bgi + Bg2f rj 
2-Btot 4 

{Bgi — Bg2){Bgl + B g2) 

2Btnt 



{P±rA 



(92) 



For production at zero rapidity, for which xi — X2 = x, 
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FIG. 14: Transverse momentum dependence of the cross section for double-gap exclusive diffraction Q. The plots show the 
modulus squared of the amplitude, iTdiffl^, with k = 1, as a function of P2±i for three fixed values of (chosen to point 
in x-direction, p'ly = 0). The upper row of plots show the contours of constant values of iTdiffj^ on a logarithmic scale as 
functions of and (in units of GeV). The lower row of plots shows the profile along the p'2y = axis. Shown are the 
results calculated with the exponential parametrizations of the two-gluon formfactor and the pp elastic amplitude, Eq. 1881 . 
The diffraction pattern in P2± evolves from a rotationally symmetric one for p'ij_ = to a two-centered one for large jp'ixl- 
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FIG. 15: Angular dependence of the cross section for double- 
gap diffractive production of a 0^ system. The plot shows 
the modulus squared of the amplitude, iTdiffp, Eq. with 
K = 1. The numbers above the curves indicate the values of 
|p'i±l — Ip'ixI- The kinematics corresponds to Higgs boson 
production at the LHC. 



and Bgi = Bg2 — Bg, this simplifies to 




(93) 



where now -Btot — B + 2Bg. In this case the amplitude 
does not depend on the variable (P^r^) = {p^^—p^^)/^, 
which is an obvious consequence of its symmetry with re- 
spect to the interchange of p^^^ and P2±- Furthermore, 
one sees that the dependences of Eq. on P]_ and 
are very different. The dependence on is monotonic 
and governed by the parameter Bg alone; it essentially 
probes the square of the two-gluon formfactors of the 
colliding protons. The dependence on r^, however, is 
governed by both Bg and B, and exhibits the diffrac- 
tive zero (I9UII : it reflects the interplay of hard and soft 
interactions. 

The qualitative differences between the P± and r± de- 
pendence of the amplitude are not specific to the expo- 
nential parametrization of the two-gluon formfactor, and 
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FIG. 16: Left: Dependence of the cross section for double-gap diffractive 0^ production on the square of the difference of 
the proton transverse momenta, r\, for P± = (i.e., Pi± = — P2X = 2rx). The kinematics corresponds to production of 
a system with M = 140 GeV at the LHC energy at zero rapidity (xi = X2 ^ 10"'^). The plot shows the modulus squared 
of the amplitude, jTdiffp, Eq. with k = 1, obtained with the exponential and dipole parametrizations of the two-gluon 

formfactor. For the exponential parametrization, the position of the diffractive zero is given by Eq. 19UII . Right: Dependence 
on the square of the sum of the transverse momenta, Pj, for r± — {i.e., — pjx ~ P±)- Also shown is Fg{t) for the two 
parametrizations, which would be the t-dependence of the amplitude without soft spectator interactions. 



can be used to test the diffractive reaction mechanism 
and extract information about the two-gluon formfac- 
tor. Figure [TBI (left plot) shows the squared moduhis of 
the ampHtude (for k — 1) as a function of r^, for P± = 0, 
in the kinematics corresponding to dijet production with 
X — 10~^ at the LHC. One sees that the dependence 
in the forward peak (near =0), and even the posi- 
tion of the diffractive dip, are not very different for the 
two parametrizations. Experimental observation of this 
structure would thus constitute a stringent test of the 
diffractive reaction mechanism. 

Figure E| (right plot) shows the dependence of the 
cross section for double-gap diffractive 0'*' production 
on Pj, for r_L = 0. Also shown is the fourth power of 
the two-gluon formfactor, Fg{t) for the two parametriza- 
tions, which would be the i-dependence of the amplitude 
without soft spectator interactions. For the exponential 
parametrization the t-dependence of the full diffractive 
amplitude is identical to that of Fg(t), cf. Eqs. and 
one sees that the two dependences are similar also 
for the dipole parametrization. The different normaliza- 
tion of the two sets of curves reflects the RGS probability, 
cf. Sec. ESI 

The position of the diffractive zero in the r'j_~ 
dependence (see Fig. ^1 left plot) is correlated with the 
slope of the monotonic P^-dependence of the diffractive 
cross section; both essentially reflect the two-gluon form- 
factor of the colliding protons. A sensible strategy for the 
analysis of double-gap diffractive dijet production would 
be to first establish the existence of the diffractive zero 
in at P± = 0, and then extract Bg from the P]_- 
dependence of the diffractive cross section at r_L = 0, 
where the cross section is maximal. In the next step. 



one could change the dijet mass (i.e., the momentum 
fractions xi = X2 = x in the two-gluon formfactor) and 
verify whether both the position of the zero and the P^- 
slope change proportionately, and whether the rate of 
change with In x is consistent with the value of a'g at the 
relevant scale, cf. Eq. H27|l. 

The a;i^2-dependence of the two-gluon formfactor in 
the protons can be probed more directly by extending 
the measurements of diffractive production to non-zero 
rapidity, j/ ^ 0, corresponding to different momentum 
fractions of the annihilating gluons in the two protons, 

= xe^y, X = (94) 

In this case the two-gluon formfactors of the two pro- 
tons are different, and the diffractive cross section is 
no longer invariant under exchange of the final proton 
transverse momenta, p'i_^_ P2-L- One sees that the 
expression for the amplitude obtained with the exponen- 
tial parametrization, Eq. (|93|l . acquires a dependence on 
(P_Lr_L) = ~p£.)/2, which is controlled by the dif- 
ference of the slopes, Bgi — Bg2. With the x-dependence 
of the individual slopes given by 

Bgi = Bg{xi) = Bg{x) + 2a'gln^, (95) 
Bg2 = Bg{x2) - Bg{x) + 2a'g\n^, (96) 

their difference is directly proportional to the rapidity 
[cf. Eq. mi 

Bg,-Bg2 = 2a'gy. (97) 
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FIG. 17: The cross section asymmetry l|98|l as a function of ra- 
pidity, y. Shown are the results obtained with the exponential 
and the dipole parametrizations of the two-gluon formfactor, 
for two representative values of a'g. In this calculation only 
the j/-dependence of the cross section arising from the con- 
volution integral is taken into account, not the y-dependence 
arising from the overall normalization, k. 

For small rapidities, \y\ ^ 1, and because of the relatively 
small value of a'g this difference is substantially smaller 
than the central value, Bg{x), and the dependence of the 
amplitude on {P±r±) can be treated in first-order ap- 
proximation. This implies that the cross section depends 
practically linearly on (Pirj_) for \y\ ^ 1. 

A convenient observable to measure is the asymmetry 
of the cross section with respect to r±_ —r±_ and Pi 
P± {i.e., P2±) fixed x and y (i.e., fixed 

xi ^ 2:2), 

0-diff(p'i_L'P2_L) + (^disiP2±,P'i±) ' 

alternatively, one could exchange the rapidities and leave 
the transverse momenta the same. This asymmetry is 
odd in y and vanishes linearly for y — > 0. For \y\ ^ 1, it is 



practically is proportional to {P±r±) = — p'2±)/2. 
When calculating the asymmetry at finite y, we have 
to take into account that, in general, also the overall 
normalization of the cross section, k{s, ^1,^,2), changes 
with y, because of the change of arguments in the t — 
gluon GPDs of the protons, see Eq. Ht)2|l . However, this 
change relative to the value at y = is of second order 
in y (the changes in the arguments of the gluon densities 
cancel each other to first order) and can be neglected for 
small y. This implies that for y <^ 1. the asymmetry is 
of the form 

A ^ Cya'giP^r^), (99) 

where the constant, C, is calculable solely from the con- 
volution integral of the two-gluon formfactor and the pp 
elastic amplitude, and does not contain information on 
the gluon densities. For finite ?/, the asymmetry is still 
proportional to a'g{P±r±), but the coefficient is a more 
complicated function of y, which depends also on the 
gluon densities in the colliding protons. Figure [T7I shows 
the theoretical estimate of the asymmetry as obtained 
with the exponential and the dipole parametrizations of 
the two-gluon formfactor, as a function of the rapidity, 
y, for two representative values of a'g. In this calcula- 
tion, for simplicity, we have taken into account only the 
y-dependence of the cross section arising from the RGS 
integral, not the y-dependence arising from the overall 
normalization, k. The latter is oc at small y but may 
be numerically important at y 1; the curves in this 
region are shown for illustrative purposes only. 

In the above discussions we have considered production 
of a O"*" system. The cross section for production of a 0^ 
state is significantly suppressed compared to 0"'". This 
is because in the hard scattering process only one gluon 
polarization state gives a large contribution in the LO 
approximation '44], and from one gluon polarization it 
is impossible to build a parity-conserving amplitude for 
the production of a 0^ state. 

Our discussion of the transverse momentum depen- 
dence in this Section is based on the approximation of 
independent hard and soft interactions, Eqs. H66H70|I . 
One should expect that the inclusion of correlations be- 
tween hard and soft interactions, as described in Sec- 
tional would modify not only the RGS probability but 
also the transverse momentum dependence of double- 
gap exclusive diffraction. This interesting question will 
be addressed elsewhere. 



VIII. DISCUSSION AND OUTLOOK 

In this paper we have outlined an approach to RGS 
in double-gap exclusive diffractive processes in pp scat- 
tering based on the idea that hard and soft interac- 
tions are approximately independent because they pro- 
ceed over widely different time-and distance scales. We 
have shown that this idea can be practically implemented 
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in the framework of Gribov's parton picture of high- 
energy scattering, and gives rise to a conccptuaUy clear 
and quantitative description of RGS. 

In the independent interaction approximation, where 
correlations between hard and soft interactions are ne- 
glected, the RGS probability can be expressed in a 
model-independent fashion in terms in two phenomeno- 
logical ingredients — the gluon GPD in the proton, and 
the pp elastic scattering amplitude. At this level we re- 
cover a simple geometric picture of the interplay of hard 
and soft interactions in the impact parameter represen- 
tation. The fact that the pp elastic amplitude at high 
energies approaches the BDL at TeV energies suppresses 
small impact parameters and ensures dominance of pe- 
ripheral scattering in double-gap diffraction. This is cru- 
cial both for justifying the approximations made in our 
derivation, and for determining the numerical value of 
the RGS probability. Our numerical results for the RGS 
probability are somewhat lower than those obtained pre- 
viously within the eikonalized Pomeron model of soft in- 
teractions; the agreement of Eq. (|83|l with the best esti- 
mate quoted in Ref. ^ is accidental and due to the fact 
that these authors assume a larger radius of the trans- 
verse spatial distribution of gluons with x 10~^. 

Our numerical prediction for the RGS probability in 
double-gap exclusive diffractive Higgs boson production 
at the LHC {mn = 100 - 200 GeV, = 14 TeV) in the 
independent interaction approximation is S'^ « 0.03. Ac- 
cording to the detailed calculations of Ref. '3| , this would 
put the estimated cross section for Higgs production in 
the minimal supersymmetric model (MSSM) with detec- 
tion through the bb mode at Bt{H bb) adis{0~^) ~ 
few X 10 fb, which should be accessible at the LHC; see 
Ref. i for details. 

Our partonic approach to RGS also allows us to discuss 
the effect of correlations between hard and soft interac- 
tions on the RGS probability. Such correlations lower 
the RGS probability compared to the independent in- 
teraction approximation. In the case of long-distance 
correlations due to the proton's pion cloud, we estimated 
this effect to be of the order of ~ 10%. Potentially more 
important are short-distance correlations, e.g. those sug- 
gested by the phenomenological concept of a "constituent 
quark" structure of the proton) . Such correlations can be 
regarded as a change of the effective size of the diffrac- 
tively scattering systems, and could reduce the predic- 
tions for the RGS probability by a substantial factor. 
In view of its importance for the Higgs boson search at 
the LHC this problem clearly requires further theoretical 
study. It can in principle also be addressed experimen- 
tally, by "measuring" the RGS probability in double- 
gap processes for which the hard scattering amplitude is 
known, such as dijet production. 

An important ingredient in our description of RGS is 
the profile function of the complex pp elastic scatter- 
ing amplitude. This underscores the importance of the 
planned measurements of pp elastic scattering and total 
cross sections in the TOTEM experiment at the LHC 



[s^ . as well as at RHIC [s^l- In addition to providing 
input for more accurate estimates of the RGS probabil- 
ity in diffraction, such measurements would allow us to 
further explore the fascinating new regime of the BDL in 
high-energy hadron-hadron scattering. 

Measurements of the transverse momentum depen- 
dence of double-gap exclusive diffractive processes with 
large cross section (dijet production) would allow one to 
perform detailed studies of the diffractive reaction mech- 
anism. Following the strategy outlined in Section I VIII 
once the reaction mechanism has been established, one 
could even use such processes to extract information 
about the transverse spatial distribution of gluons in the 
colliding protons, including its change with x. Such stud- 
ies would complement the information on the two-gluon 
formfactor obtained from vector meson production at 
HERA or a future electron-ion collider (EIC). Eventu- 
ally, using QCD evolution as well as models of nucleon 
structure, these data on the transverse spatial distribu- 
tion of gluons could also be correlated with the planned 
measurements of quark GPDs in hard exclusive processes 
in ep scattering in fixed-target experiments (HERMES, 
JLab 12 GeV, COMPASS). One of the advantages of 
Gribov's parton picture of hard and soft interactions 
is precisely that it unifies the description of hadron- 
hadron and electron/photon-hadron scattering at high 
energies. Other ways to probe GPDs in pp scattering 
with hard processes (non-diffractive) have been described 
in Ref. Tq. 

We would like to comment on some of the experimen- 
tal aspects of measurements of the transverse momentum 
dependence of double-gap exclusive diffraction with the 
proposed forward detectors at the LHC. Such measure- 
ments require good energy resolution in order to guaran- 
tee exclusivity and determine the mass of the diffractively 
produced system, as well as sufficient transverse momen- 
tum resolution to map the P2_L distributions. An im- 
portant experimental problem is that the intrinsic trans- 
verse momentum distribution in the beams at the interac- 
tion point (IP) puts a lower bound on the transverse mo- 
mentum transfers that can be resolved. This distribution 
is determined by the beam optics, and thus closely corre- 
lated with the luminosity. The proposed 420 m forward 
detectors for the CMS and ATLAS experiments [H HI 
can tag protons in the range 0.002 < ^ < 0.015; in the 
TOTEM experiment at CMS with detectors at 200 m the 
range will be extended to ^ < 0.1 [s^- Both detectors 
can obtain a longitudinal and transverse momentum res- 
olution comparable to the beam distributions. At a lumi- 
nosity of lO'^'^ cm~^s~^ with (3* = 0.5 m, and a one-sigma 
normalized emittance e = 3.75 x 10^^ tt m, the one-sigma 
angular spread of the beams at the IP is 8 fj,r, correspond- 
ing to a transverse momentum spread of 56 MeV/c. This 
sets the scale for experimental smearing of the transverse 
momentum distributions. The TOTEM experiment also 
envisages running at substantially larger (3* values (18, 
90 and 1540 m). These values will reduce the transverse 
angular spread of the beams at the IP by yJe/{iTj3*), but 
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with a concomitant reduction in luminosity. Given that 
the typical scale of the transverse momentum distribu- 
tions is TTig w 1 GeV, it seems feasible to make detailed 
and precise measurements of the transverse momentum 
distributions with the LHC420 and TOTEM detectors 
even when running in high-luminosity mode. 
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[63] In the present context it is convenient to define the gluon 
GPD as corresponding to the gluon momentum density, 
xG(x), not the number density, G{x), in the limit of zero 
momentum transfer. 

[64] The Sudakov formfactor associated with the ggH vertex 
accounts for the absence of gluon bremsstrahlung with 
transverse momenta |fcj_| < |fe_L,rad| < mn l^,, radiated 
from the annihilating gluons. The requirement of absence 
of radiation from the screening gluon results in an addi- 
tional shift of the k± distribution towards larger values. 
However the amplitude for radiation from this line does 
not contain factors ois \n{k\ / m}fj\&nd is beyond the ac- 
curacy of the calculation of Ref. j35| . 
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integral of squared spatial distribution of gluons) is differ- 
ent from the distribution Pi{b) (square of overlap integral 
of spatial distribution of gluons), which was introduced 
in Ref. [3 describe double dijet production in inclu- 
sive high-energy pp/pp scattering. The two distributions 
coincide only in the case of an exponential parametriza- 
tion of the two-gluon formfactor, which corresponds to 
a Gaussian dependence of the formfactors on transverse 
momenta and coordinates; see Eqs. (EU and ((THJ. In gen- 
eral, it is not correct to replace Phard(b) by Piib), as was 
done in Refs. ESIIl. 

[66] The value 5*^ = 0.026 quoted in Ref. is obtained when 
taking into account corrections to the hard scattering am- 
plitude resulting from the proton transverse momenta, 
which are not included in our approach. The difference 
between the two values is immaterial for the present dis- 
cussion. 
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not result from the presence of the fluctuations of the 
cross section strength. As a result, the cross section of 
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